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Abstract 

For an ideal / C M.[x] given by a set of generators, a new semidefi- 
nite characterization of its real radical /(Mr(/)) is presented, provided 
it is zero-dimensional (even if / is not). Moreover we propose an al- 
gorithm using numerical linear algebra and semidefinite optimization 
techniques, to compute all (finitely many) points of the real variety 
Vr(/) as well as a set of generators of the real radical ideal. The latter 
is obtained in the form of a border or Grobner basis. The algorithm is 
based on moment relaxations and, in contrast to other existing meth- 
ods, it exploits the real algebraic nature of the problem right from the 
beginning and avoids the computation of complex components. 
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1 Introduction 

Algebraic computation over the reals is a highly relevant topic with many 
practical applications and, in particular, for finding real solutions to a sys- 
tem of polynomial equations. Throughout the paper, K[x] = K[xi, . . . ,Xn] 
denotes the ring of polynomials in n variables over the field IfC = M or 
C. For an ideal I C K[x], Vc{I) := {x G | f{x) = V/ G /} and 
V]r(/) := Vc{I) n denote, respectively, the complex and real varieties of 
/ and, for V C C", I{V) := {/ G K[x] \ f{v) = e V} is the vanishing 
ideal of the set V. The ideal I{Vc{I)) coincides with the radical ideal \/7 
of / by the Nullstellensatz and I(Vk{I)) coincides with the real radical ideal 
by the Real Nullstellensatz (see Section 2.1 for details). The problem of 
finding the radical ideal /(Vc(/)) seems to be much better understood than 
that of finding the real radical ideal /(Vr(/)); see below for a brief recap 
on existing literature. In this paper, we provide a new characterization of 
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the real radical ideal I(Vk(/)) of an ideal / C W[x], assuming / is given by 
generators hi, ... , hm G M.[x\ and Vm.{I) is finite (while Vc(/) needs not be 
finite). In addition, from this characterization, we also define a numerical 
algorithm based on semidefinite programming to compute the points of the 
(finite) variety V^{I) as well as a set of generators of the real radical ideal 
I(y^{I)). More generally, our results extend to the case of the so-called 
^-radical ideal I{V^{I) fl S) where S C R" is defined by finitely many poly- 
nomial inequalities, assuming that MR(-f) n S is finite. It turns out that 
a similar algorithm also works for computing Vc(/) and the radical ideal 
/(Vc(I)) (assuming now Vc{I) finite) although very good methods already 
exist for this latter case. In the remainder of the Introduction, after recall- 
ing some motivation and related literature on the problem of finding the 
(real) radical ideal, we sketch the main ingredients of our method. We al- 
ready introduce some definitions but refer to Sections 2 and 3 for additional 
definitions about polynomials and moment matrices. 

Motivation. The main motivation of this work is to provide a characteriza- 
tion as well as an algorithm for finding the real variety and the real radical of 
an ideal I C R[x] that takes into account the specific real algebraic geomet- 
ric nature of the problem. Indeed, to the best of our knowledge, most basic 
methods for computing the real variety Vk(/) first compute the complex 
variety Vc (/) ; for this they require as basic ingredients a Grobner basis of / 
and a linear basis of the vector space M[x]// and thus they work under the 
assumption that V(c{I) is finite. Even if Vc(/) is finite but has many more 
complex elements than real ones, this may produce a large computational 
overhead. This is particularly important as the numbers of complex and real 
solutions may differ significantly as supported by the fewnomial theory of 
Khovanski [22]; see also the discussion in Bihan et al. [4], [5]. In other words, 
this problem of real algebraic geometry is solved via algebraic methods that 
do not take into account right from the beginning the real algebraic aspect 
of the problem. In contrast, our characterization and our algorithm do not 
need knowledge of a Grobner basis of I and are real algebraic in nature, as 
we never compute any complex zero. 

Related literature. There is a large literature on the problem of finding the 
radical ideal y/l of an ideal /; see, e.g., [3], [8], [17], [19], [23], [24]. For the 
general (positive-dimensional) case, Krick and Logar [23] propose an efficient 
algorithm based on splitting and reduction to the zero-dimensional case, 
which is implemented e.g. in the computer algebra package Singular [18]. 
In the zero-dimensional case the problem is considered to be well-solved, e.g., 
via the following method of Seidenberg [37]: ^/l = (/ U {qi, . . . , qn}), the 
ideal generated by / and the qis, where qi is the square-free part of the monic 
generator pi of /nK[xj]. Finding pi is easy once a linear basis B of ]K[x]// is 

known. Namely, find the smallest integer ki for which {l.Xj. xj .t*^'} is 

linearly dependent in ]IC[a;]//; then this smallest linear dependence gives the 
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polynomial pi. Next, the polynomial qi can be found taking derivatives and 
gcd-computations as qi = So finding /(Vc(/)) is easy if we have a 

basis of ]K[a:]//. A classical method for finding such a basis B is to compute 
a Grobner basis of I and the corresponding set B of standard monomials. 
The results in the present paper show that one may alternatively find such 
a basis B from a suitable moment matrix. 

On the other hand, the problem of computing the real radical ideal is 
considered to be much more difficult. For instance, in their paper which 
is one of the first classical references on this problem, Becker and Neuhaus 
[1, p. 7] write that the computation of T-real parts (thus, the real radical 
ideal) is much more difficult (than that of the ordinary radical). They give 
an algorithm for based on finding the minimal real prime ideals Pi such 
that ^/l = rijPj. Among other advanced algebraic manipulations, their 
algorithm makes intensively use of (ordinary) radical computations. For 
other works along similar lines see, e.g., [2], [9]. 

Finally, let us mention that excellent algorithms and software packages 
exist for computing the complex variety Vc (/) of a zero-dimensional ideal /, 
e.g., by Verschelde [44], Rouillier [36]; see also related work by Mourrain et 
al. [32] and e.g. the monograph [16]. For instance, Verschelde [44] proposes 
symbolic-numeric algorithms via homotopy continuation methods (cf. also 
[38]) whereas Rouillier [36] solves a zero-dimensional system of polynomials 
by giving a rational univariate representation (RUR) for its solutions, of the 
form f{t) = 0, xi = . . ., x„ = where f,g,gi, ■■■,gn G ^[t] are 

univariate polynomials. The computation of the RUR relics in an essential 
way on the multiplication matrices in the quotient algebra ]K[a;]// which 
thus requires the knowledge of a corresponding linear basis. 

Our contribution. Given an ideal I C K[x] (K = R, C) defined by a set of 
generators and satisfying |T^(7)| < oo, we provide a method for computing 

Vk{I) as well as a border basis and a Grobner basis of the ideal I{V^{I)). 
Our approach is based on a semidefinite programming characterization of 
/(Vk(/)) with the following distinguished feature. Remarkably, all infor- 
mation needed to compute the above objects is contained in the so-called 
moment matrix (whose entries depend on the polynomials generating the 
ideal /) and the geometry behind it when this matrix is required to be pos- 
itive semidefinite with maximum rank. The latter property is achieved by 
standard semidefinite programming algorithms. For the task of computing 
the real roots and the real radical ideal I{Vm.{I)), the method is real alge- 
braic in nature, as we do not compute (implicitly or explicitly) any complex 
element of Vc(/). 

Lasserre [25] already recognized that moment matrices can be used for 
approximating the minimum of a polynomial over a basic closed semi-algebraic 

set and sometimes extracting global minimizers (cf. [21]). The present paper 
builds on this approach and shows how it can be applied to finding the real 
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radical of a zero-dimensional ideal. Moreover there are links between mo- 
ment matrices and the Hermite quadratic forms used in [33] for computing 
the number of real roots, that were pointed out in [28]. 

Our approach with its specificity is best illustrated on the task of com- 
puting the real radical ideal I{Vr{I)) that we now briefly describe. 

Given a sequence y = {ya)a€N" € M^", consider the moment matrix 

(later we will also introduce complex moment matrices M'^{y), M^'^(y)). 
One may think that y and M^{y) are indexed by the set T„ := {x" | a G N"} 
of monomials. Given a polynomial h € R[x], set vec{h) := {ha)am" and 
define the new sequence hy := M^{y)vcc{h) G M^". By abuse of language let 
us say that h lies in the kernel of M^{y) when vec(/i) does, which enables us 
to view KerM''*(y) as a subset of R[x]. The following property of moment 
matrices plays a central role in our approach; it is based on ideas from 
[12], [13], [27] and will be proved at the end of Section 3. Let I = {hi, . . . , hm) 
be an ideal generated by /ii, . . . , hm G M[x]. 

Proposition 1.1. Assume that Vk(/) is finite. If 

M^(y) y 0, M'^ihjy) = (i = 1, . . . , m) (1.1) 

then the kernel of M^(y) is a real radical ideal, rankM'^(y) < |T^(/)| and 
/(Vr(/)) C KevM*{y), with equality if and only if M^{y) has maximum 
rank, equal to |T^(/)|. 

(In (1.1) the notation 0" stands for positive scmidcfinitc.) This semidefi- 
nite characterization leads directly to an algorithm for computing /(Vm(/)), 
by considering truncated moment matrices in place of the full (infinite) mo- 
ment matrix M^{y). Namely, given an integer t, let Mf-(y) denote the 
principal submatrix of M*{y) whose rows and columns are indexed by the 
set Tn,t '■= {x" I a G N" with |a| := Yli (^i ^ t} and set 

dj := [deg{hj)/2], d := max dj. (1.2) 

j=l,...,m 

Fix t > d and assume M^{y) is a maximum rank matrix satisfying 

Mf (y) y 0, Ml^. {hjy) = (j = 1, . . . , m). (1.3) 

We will show that if, moreover, 

rankMf (t/) = rankMf_rf(y) (1.4) 

for some d < s < t, then /(Vk(/)) coincides with the ideal generated by 
KeTMf-{y). The same conclusion holds if 

rankiVf,^(y) = iankMf_^{y) (1.5) 
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for some 2d < s < t. Moreover, from the semidefinite characterizations 
(1.3)-(1.5), the fohowing algebraic objects can be obtained directly from the 

matrix Mf (y): 

(i) Let B C Tn^s be a set indexing a maximum nonsingular principal sub- 

matrix of Mf{y). Then ;B is a linear basis of the quotient vector space 
R[x]/I{Vm.{I)) (see Section 3.3). 

(ii) We can compute directly from M^{y) the matrix of any multiplica- 

tion operator in R[x]/I{Vm{I)) with respect to the basis B, and thus 
compute Vm,{I) (using the eigenvalue method, see Section 2.2). 

(iii) When the set B (as in (i)) is an order ideal (i.e., is stable under divi- 
sion), the matrices of the multiplication operators by xi, . . . ,Xn give 
directly a border basis of the ideal I(V^(/)) (see Section 2.5). 

(iv) Given a graded lexicographic monomial ordering, we can find a set 
B (as in (i)) which is precisely the set of standard monomials; the 
associated reduced Grobner basis of /(Vk(I)) can then be recovered, 
since it is contained in the border basis. In fact our method also applies 
to an arbitrary monomial ordering (see Section 4.4.5). 

(v) Finally the method can also detect whether the real variety T^(/) is 

empty. Indeed, Vk(/) =0 if and only if, for some integer t, the system 
(1.3) admits no solution y with yo / 0. (See Remark 4.8.) 

Further discussion. An independence oracle in R[x]/I{V]s_{I)) is needed 
for our algorithm in (iv) above. The following property is a crucial ingredi- 
ent. Assume that one of the conditions (1.4) or (1.5) holds and consider a 
set T C Tn^s- Then, T is linearly independent in R[x]//(Vr(I)) if and only 
if T indexes a linearly independent set of columns of M^{y). In view of (iv), 
a Grobner basis can easily be derived afterwards in contrast with classical 
methods which compute the set of standard monomials from the Grobner 
basis. 

Realizing the above tasks relies only on numerical linear algebraic oper- 
ations on Mf^{y) like evaluating the rank of certain principal submatrices. 
Finding a matrix satisfying (1.3) is an instance of semidefinite programming. 
Moreover, it is a property of most interior-point algorithms for semidefinite 
programming that they do find such a matrix having maximum rank (see 
Section 4.4.1 for details). 

The method is iterative. Namely, if the maximum rank matrix satisfying 
(1.3) does not satisfy (1.4) or (1.5), then we iterate with t -|- 1 in place of t. 
The method eventually terminates since we will show that (1.4) holds for t 
large enough. 

The following two small examples illustrate how positive semidcfiniteness 
of the matrix Mf^{y) allows the elimination of all complex (nonreal) roots. 
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whose number can be much larger than the number of real roots or even 
infinite. 

Example 1.2 Let / C R[x] be generated by hi = Xi{xf + 1) (i = 1, ..,n). 
Then, Vr{I) = {0}, |Vc(/)| = 3", di = 2 for all i. Assume y satisfies (1.3) 
for order t = 3. Then Mf{hiy) = implies y4e. = -?/2ei and Mf{y) ^ 
implies y2ei-,yiei > which in turn implies = for all a 7^ with |a| < 5. 
(Throughout, ei,...,en denote the standard unit vectors in M".) Hence 
rankM|^(y) = rankM^(y) = 1; that is, (1.4) holds for s = 2. In fact, 
KerM]^(y) is spanned by xi, . . . the generators of I(Vr(I)). One may 
argue that the ideal I is already described by a Grobner basis. But the same 
conclusion also holds under the change of variables x = Ay with A being a 
nonsingular matrix, in which case other methods would require a Grobner 
basis computation. 

Example 1.3 Let I C R[xi, a;2] be generated hy h = x\+x\. Then Mr(/) = 
{0} and Vc(/) = {(xi,X2) | xi = =bia;2} is infinite. Then M^{hy) = gives 
2/2ei + 2/262 = which, together with Mf{y) y 0, implies = for a 7^ 0. 
Hence the maximum rank of M^{y) is equal to 1 and KerMf (y) is spanned 
by xi,aj2, the generators of /(Vk(/)). 

The method sometimes (partially) applies even if none of the rank con- 
ditions (1.4), (1.5) holds, namely when Mf\y) contains sufficient informa- 
tion for the construction of the (formal) multiplication matrices. More pre- 
cisely, let Mf-{y) be a maximum rank matrix satisfying (1.3), let B C T„^t 
index a maximum nonsingular principal submatrix of M^{y), set dB := 
(U^^iXiB) \ B, and assume that the two principal submatrices of Mf-{y) 
indexed by B and by ;S U dB have the same rank. Then, by the results of 
Kehrein, Kreuzer and Robbiano [16, Ch. 4], we can construct the formal 
multiplication matrices and, if they commute pairwise, asetW^ Vm{I) can 
be computed. By checking whether the points of W satisfy all the equations 
hj = 0, we can eliminate the points in W \ Vm.{I). It turns out that, for 
most examples we have tested, W = Vm.{I) holds and we are again able to 
find I{Vm.{I) together with a border basis generating this ideal. 

Finally, the method also applies to the task of finding the radical ideal 
I{Vc{I)) of a zero-dimensional ideal / C C[x]. For this, instead of using 
the matrix Mf-{y) where y is a real sequence indexed by Tn,t, we have to 
use a matrix Mf^iy) where the argument is a complex sequence indexed by 
^2n.t (sec Section 3.1 for details). Similar results hold as in the real case. 
Namely, under certain rank conditions, the ideal /(Vc(/)) can be obtained 
as the ideal generated by the kernel of a maximum rank complex moment 
matrix (see Section 4.3 for details). However, a drawback in the complex 
case is that one must in general handle matrices of larger order which leads 
to larger semidefinite programs, thus more difficult to solve. However, so 
far we do not claim that our method can compete with existing methods for 
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finding the complex variety Vc{I) as e.g. [36], or [44], especially in view of 
the present status of SDP solvers (that we use as a black box), still in their 
infancy. 

Contents of the paper. Section 2 contains preliminaries about ideals of 
polynomials, in particular, about the quotient ring ]K[x]//, multiplication 
matrices, Grobner bases and border bases. We also indicate in Section 2.4 
an algorithm for finding the set of standard monomials from an independence 
oracle in ]K[x]//. Section 3 contains preliminaries about moment matrices, 
in particular, results relating (real) radical ideals and kernels of positive 
semidefinite moment matrices. In Section 4, we prove the main results about 
the semidefinite characterization of the variety Vk{I) and the associated 
radical ideal I{Vk{I)). Section 4.4 gives the details and implementation 
of an algorithm based on the semidefinite characterization, and Section 5 
contains several examples illustrating its behaviour. 

2 Preliminaries on Polynomial Ideals 
2.1 Polynomial ideals and veirieties 

Throughout, IfC = M or C, and K[x] := lK[xi, . . . , x„] denotes the ring of 
multivariate polynomials in n variables over the field K. For an integer 
t > 0, K[x]t denotes the set of polynomials of degree at most t. For a scalar 
a G C, a denote its complex conjugate and, for a vector u e (resp., a 
matrix A), u* (resp.. A*) denotes its conjugate transpose. Following e.g. 
[10], denotes the monomial x"^ • • • x"" (for a G N") and cx" is called a 
term (for c G K.). Let T„ := {x" | a G N"} denote the set of monomials 
and set NJ^ := {a G N" | \a\ := Y.7=i «i < *}, T„,t := {x" | a G N?} for 
i G N. Following [16, Ch. 4], a set B C T„ is called an order ideal if B is 
stable under division, i.e., for all a,b G T„, b G B, a\b implies a G B. Given 
an ideal"^ / C K[x], let 

Veil) := {x G I fix) = V/ G /}, Fm(/) := Vc{I) n 

denote its complex and real varieties, respectively. For a set F C K", define 
the ideal 

I{V) := {/ G K[x] I f{v) = 0\/ve V}. 
Given an ideal / C K[x], one can define the ideals /(Vc(/)) and 

V7 := {/ G K[x] If"" el for some m G N \ {0}} 

''An ideal / C K[x] is an additive subgroup of K[x] such that fg € I whenever f € I 
and g € K[x]. 
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and, when I C M[a;], one can define the ideals I{y^{I)) and 

Vl ■.= {p(^ M[x] I p^™ + ^ ?| G / for some G M[x], rn G N \ {0}}. 

j 

Obviously, 

I^VlC I{Vc{I)), ICVIC I{Vm{I)). 

The ideal J is said to be radical (resp., real radical) if I = s/l (resp., I = ^fl). 
Obviously, / C I(Vc{I)) C I{Vr{I)). Hence, if / C R[x] is real radical, then 
/ is radical and moreover, Vc{I) = Vk.{I) ^ ]R"if < oo. The following 

lemma gives a useful characterization for (real) radical ideals. 

Lemma 2.1. An ideal I C IfC[a:] is radical if and only if 

VpeK[x], p'^el^pel (2.1) 

and / C R[x] is real radical if and only if 

Vpi,...,pfc gM[x], pl + ...+plel^pi,...,Pkel. (2.2) 

Proof. If / is radical, then (2.1) obviously holds. Conversely assume 
that (2.1) holds; we show that p'^ e I =^ p G / by induction on m > 2. If 
p"^ e I then G / and thus, by (2.1), pT^/^l £ / which, by the induction 
assumption, implies p £ I. The proof for the real radical case is along the 
same lines and thus omitted. □ 

Theorem 2.2. (i) Hilbert's NuUstellensatz (see, e.g., [10, U-1]) Vl = 
I{Vc{I)). 

(ii) Real NuUstellensatz (see, e.g., [7, %4-l]) = I{yM.{I)) for an ideal 
I C M.[x\. 

For a polynomial p G x 1— > p{x) = ^cePa^^^ vcc{p) := (f'a)aeN" 
denote the vector of its coefficients. We also let vec(p) denote the vector 
{Pa)a€f^'^ for any t > deg(p), as Pa = whenever |a| > deg(p). 

Finally, given A C K[x], let (A) := {J2i=iUiPi \ Pi e A,Ui e K[x]} 
denote the ideal generated by A. 

2.2 The algebra ]K[a;]// and multiplication matrices 

Consider the quotient space K[x]/I, whose elements are the cosets [/] := 
f + I = {f + q\ q G /} for / G K[x]. K[x]/I is a K- vector space with 
addition [/] + [g] := [/ + g] and scalar multiplication A[/] := [A/], and 
an algebra with multiplication [f][g] := [fg], for A G K, /, 5 G M.[x]. In 
particular, for h G K[x], the multiplication operator 

ruh : K[x]/I — > K{x]/I 
if] ^ [hf] 
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is well defined. The following well known result relates the cardinality of 
Vc{I) and the dimension of the vector space ]K[x]//. See e.g. [10], [40] for a 
detailed treatment of the quotient algebra ]K[.x]//. 

Theorem 2.3. For an ideal I in K[x], \Vc{I)\ < oo <^=^ dimIK[a:]// < oo. 
Moreover, \Vc{I)\ < dim K.[a;]//, with equality if and only if I is radical. 

Assume |Vc(/)| < oo and set N := dimK[a;]/I > |Vc(-^)|- Consider a 
set B := {bi, . . . , bi\f} K[x] for which the cosets [fei], . . . , [^at] are pairwise 
distinct and {[bi], . . . , [^at]} is a basis ofK[x]/r, by abuse of language we also 
say that B itself is a basis of ]IC[x]//. Then every / G IfC[x] can be written 
in a unique way as / = Xli^i '^i^i + where Cj G K, jo G /; the polynomial 
rcsi^{f) := Cibi is called the residue of f modulo I w.r.t. the basis 

B. In other words, the vector space Span]jj(i3) := {X^^i (hbi | Cj G K} is 
isomorphic to K[x]/I. 

Given a basis B of ]K[a;]// and h G K[x], let Xh denote the matrix of the 
multiplication operator with respect to B. That is, writing resis{hbj) = 
'^iLi dijbi, the jth column of is the vector (a^j)^^. The following well 
known result relates the points of the variety Vc{I) to the eigenvalues and 
eigenvectors of X^. See, e.g., [16, Ch. 2,3] for a detailed treatment. 

Theorem 2.4. Let h G K[x] and, for v G Vc{I), set Cb,v ■= {bi{v))fLi. The 
set {h{v) I V G Vc(-^)} is the set of eigenvalues of Xu o.'nd X^(b,v = h{v)C,B,v 
for all V G Vc(/). 

When the matrix is non-derogatory (i.e., all its eigenspaces are 1-dimensional), 
one can recover the points v ^V£{I) from the eigenvectors of X^ . If / is rad- 
ical, then = |Vc(/)| and thus Xy^ is non-derogatory whenever the values 
h{v) {v G Vc{I)) are pairwise distinct. This is achieved with high probability 
if one chooses h = Yl^=i ^i^i random scalars Oj. 

2.3 Grobner bases and standard monomials 

A classical basis of IK[x]// is the set of standard monomials with respect to 
some monomial ordering V of T„. Let us recall some definitions. (See e.g. 
[10] for details.) Fix a monomial ordering >- on T„. Write also ax" >- bx^ 
if x" >- x^ and a,b G IK \ {0}. For a nonzero polynomial / = Xla/o^"' 
its leading term LT(/) is the maximum /ax" with respect to >- for which 
/„ / 0. The leading term ideal of / is LT(I) := (LT(/) | / G /) and the set 

By := T„ \ LT(/) = {x° | LT(/) does not divide V/ G /} 

is the set of standard monomials. Obviously By is an order ideal. A finite 

set G C / is a Grobner basis of I if LT(/) = {LT{g) \ g e G); thus x" G By 
if and only if .t" is not divisible by the leading term of any polynomial 
in G. A Grobner basis always exists and it can be constructed, e.g., with 
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the algorithm of Buchberger. Call G reduced if, for all g G G, the leading 
coefficient of hT{g) is 1 and no term of g lies in {LT{g') \ g' G G \ {g})- 
Given nonzero polynomials f,hi,..., h^, the division algorithm applied to 
dividing / by /ii, . . . , hm produces polynomials ui, . . . , Um, r satisfying / = 
X^jLi Ujhj+r, no term of r is divisible by LT{hj) (j = 1, . . . , m) and LT(/) >z 
LT{ujhj). Note that deg(uj/ij) < deg(/) when the monomial ordering is a 
graded lexicographic ordering. When {hi, . . . , hm} is a Grobner basis of the 
ideal I := {hi, ... , hm), the remainder r is uniquely determined and belongs 
to Spanjj(i3^); moreover, / G i" r = 0. Therefore, the set By is a basis 
ofK[x]/I. 

For an arbitrary basis B of IC[a:]//, set := maxjgg dcg(?)). The next 
result shows that ds is minimum when B is the set of standard monomials 
for some graded lexicographic order. 

Lemma 2.5. Let I be a zero-dimensional ideal in K[x]. Let {gi, ■ ■ ■ ,gk} 
he the Grobner basis of I with respect to a graded lexicographic monomial 
ordering and let By be the corresponding set of standard monomials. For 
any basis B ofK[x]/I, we have ds^ < d^- 

Proof. Set B = {bi, bN}- Write hi = Ex«gb^ <^,ocX°' + Ylh=i '^h.gh 
where q^q, G K for i = 1, . . . , and G IK[x]. Then deg{uhgh) < deg(6j) 
(by the properties of the division algorithm as we use a graded mono- 
mial ordering). Thus, deg(J2x'^eB^ Ci^ax") < deg(6i). Let a;°° G By with 
deg(a;°°) = dsy As B, By are two bases of K[a;]/7, the matrix {ci^a) i=i,...,N 

aeBy 

is nonsingular and thus its aoth column is nonzero. Hence Ci^ao for some 
i. Hence ds^ = deg{^^c.^By Ci,aX°') < deg(6i) < ds- □ 

2.4 Finding the set of standard monomials from an indepen- 
dence oracle 

When / is a zero-dimensional ideal and ;^ is a monomial ordering on T„, we 
describe a method for finding the set By of standard monomials, assuming 
we have an oracle for checking linear independence in ]K[a;]/I. This 'greedy 
sieve' algorithm, described below in Algorithm 1, does not require knowledge 
of a Grobner basis of /. 

The next lemma shows correctness of Algorithm 1 and how to use it for 
finding the set By of standard monomials. 

Lemma 2.6. Let I be a zero-dimensional ideal, >- a monomial ordering on 

Tn, and By = T„\LT(/), the associated set of standard monomials. For an 
integer s > 1, let Bg be the set returned by the greedy sieve algorithm applied 
to il,y,s). 

(i) Bg is linearly independent in K[x]/I and satisfies By H T„_s C Bs; in 
particular, Bs = By if By C T^^s- 
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Algorithm 1 The 'greedy sieve' algorithm: 

Input: A zero-dimensional ideal / G M.[x], a monomial ordering >- on T„, 

and an integer s > 1. 
Output: A set B C T^^ linearly independent in IK[x]/I and satisfying 
B^By^] Tn,s 
1: Order the monomials in Tn,s with respect to y. 
2: Initialize B := $, L := (^1,^2, • • ■), the ordered set Tn,s- 
3: while B C L do 

4: Set t as the first element of L\B 

5: if ^ U {t} is linearly independent in then 

6: Reset B:=Bu{t} 

7: else 

8: Reset L := L \ tT„ (i.e., remove from L all multiples of t). 

9: end if 
10: end while 
11: return B = L 



(ii) IfBs = Bs+i thenBs=By. 

(iii) If>- is a graded monomial ordering, then Bg C B-^ ; therefore, Bg = By 
if\Bs\ =dimK[a;]//. 

Proof, (i) Obviously, throughout the algorithm, B is linearly independent 
in and B C L. Assume tk G {By fl T„^s) \ Bg- Consider the step 

when the algorithm examines and let B be the current set maintained 
by the algorithm. Then, B C {fi. . . . ,tk-i}, tk E L \ B and B U {tk} is 
linearly dependent in ]K[x]//. Hence there exists a polynomial f E I with 
LT(/) = tk, contradicting the assumption that tk ^ By. This shows By fl 
Tn,s ^ Bs- Moreover, if By C T^^^, then By C Bg] equality holds since 
\Bg\ < dim IK [a;] /I as Bg is linearly independent in while \By\ = 

dimIK[a;]/I. 

(ii) Assume Bg = Bg+i. Then, in view of (i). By n {Tn,g+i \ T„,s) = 0. This 
implies By C Tn^g- Indeed assume t G By has degree at least s + 1; then 
any divisor t' of t with degree s + 1 lies in By (since By is an order ideal) 
and thus t' G By Pi (T„^s_|_i \ Tn,g) = 0, a contradiction. Therefore, by (i), 

Bg=By. 

(iii) Assume is a graded monomial ordering and, say. By C Tn^d for some 
integer d. li d < s then Bg = By hy {}). If d > s + 1, then Bg C Bd (since 
all elements of l^n,d \ Tn,* come after the elements of T„^s in the ordering )^) 
and Bd = By (by (i)), implying BgCBy. □ 

Remark 2.7. Observe that, when >- is not a graded monomial degree or- 
dering, one cannot claim the inclusion Bg C By. For instance, consider the 
ideal / = (.7r^ — 1, — y + .x^ + x +1) in y] and choose as monomial order- 
ing >- the lexicographic order with y > x. Then, By = {l,x,x^} is found 
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when applying Algorithm 1 to (/, s = 2); observe that 82 = 63. However, 
the algorithm applied to {I,y,s = 1) returns the set Bi = {l,x,y}; thus 
Bi^By, while \Bi\ = 3 = dimR[x]/I. 

Alternatively, one could initialize the set L in Algorithm 1 to be the full 
ordered set T„. Then the algorithm still terminates in finitely many steps 
(beeause dimIfC[a:]// < 00) and the set B returned by the algorithm is equal 
to By (using the same argument as in Lemma 2.6 (i), one can show that 
By C B, implying By = B). 

A crucial tool for applying Algorithm 1 is having an oracle for testing 
linear independence in IK[a:]/I. In our setting the oracle will work as follows: 
Given a subset B C T„^s, B is linearly independent in K[a;]// if and only if B 
indexes a linearly independent set of columns of a suitable moment matrix 
Ms{y). This motivates why in our presentation of Algorithm 1 we explore 
the set Tn,s of monomials of degree at most s. 

2.5 Border bases and formal multiplication matrices 

We recall results about border bases following the exposition from [16, Ch. 4]. 
See also [40] for details about border bases. Given an order ideal B C T„, 
the border of B is the set 

dB := {xix^ Ix^^ eB, i = 1, . . . ,n} \ B. (2.3) 

Assume 6/0, set N := \B\, H := \dB\ and write B = {61, . . . ,6Af} and 
dB = {ci, . . . ,ch}- a set of polynomials G = {gi, . . . ^qh} is called a B- 
border prebasis if each gj is of the form 

N 

gj = Cj — ^^aijbi for some aij G K. (2-4) 

i=l 

One also says that gj is marked by the element Cj of dB. Given a polynomial 
/, the border division algorithm [16, Prop. 4.2.10] produces polynomials 
Uj,r such that / = X^^i Ujgj + r, and r G Spanjj(i3). Hence, for any ideal 
I containing G, B spans the K- vector space K[x]//. The set G C 7 is said 
to be a B-border basis of I if B is linearly independent in K[x]//, i.e., if B is 
a linear basis of K[a;]//; in that case G generates the ideal /. 

Stetter [40] advocates using border bases instead of Grobner bases since 
they do not depend on any monomial ordering. Border bases represent in 
fact an extension of the notion of Grobner bases. Indeed, the set Tn\ B 
defines a monomial ideal; the elements of the minimal set of generators of 
this monomial ideal are called the corners of B, which belong to dB. When 
B = By is the set of standard monomials for some monomial ordering, there 
exists a unique By-horder basis G of I and the reduced Grobner basis of / 
is the STibset of G consisting of the polynomials in G that are marked by the 
corners of By . 
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When G is a B-border prebasis, one can mimic the construction of the 
multiphcation matrices from the previous section in the following way. Fix 
G {1, . . . ,n}. The formal multiplication matrix is the N x N matrix 
whose zth column is defined as follows. If Xkhi G say, Xkbi = h^, then the 
zth column of is the standard unit vector e^. (with all zero entries except 
1 at the rth position). Otherwise, Xkbi € dB, say, x^hi = Cj, then the ith 
column of X^ is the vector {aij)^^ (compare with Eqn.(2.4)). We will use 
the following result (see [16, Thm. 4.3.17]). 

Theorem 2.8. Let B C.Tn be an order ideal, let G be a B-border prebasis 
with associated formal multiplication matrices Xi, . . . ,Xn, and let J := {G) 
be the ideal generated by G. Then, G is a border basis of J if and only if 

the matrices Xi, . . . ,Xn commute pairwise. In that case, B is a linear basis 
ofK[x]/J and the matrix X^ represents the multiplication operator rux^. of 
K[x]/J with respect to the basis B. 

Remark 2.9. Following Mourrain [31], call B CTn connected to 1 if 1 e B 
and any monomial in B is of the form Xi-^^Xi^ " ' ^ik with a;^^, x^^Xij, • • • , 
Xj^Xjj • • • a^ifc € B. Obviously if B is an order ideal then B is connected to 
1. As shown by Mourrain [31, Th. 3.1], the result of Theorem 2.8 remains 
valid in the more general setting where B is connected to 1 (instead of being 
an order ideal). We restrict our attention in this paper to monomial bases 
of K[x]/J that are order ideals, in particular, because we have an algorithm 
for finding such bases, as we just saw in the preceding section. It will be 
interesting to investigate the use of bases satisfying Mourrain's criterion in 
subsequent work. 

3 Preliminaries on Moment Matrices 

3.1 Moment matrices 

Given a sequence y G M^", its real m,oment matrix M^{y) is the real sym- 
metric matrix indexed by N" whose (a,/3)th entry is ya+p-, for a,P E N". 
Given a sequence y G its complex moment matrix is the matrix M'^'^iy) 

indexed by N^" whose (aa', /?/?')th entry is ya'+p^a+p' ■, for ("jQ!')' W^P') ^ 
N^". If y satisfies 

Va'a = W for (a, a') G N^", (3.1) 

then M^'^iy) is a Hcrmitian matrix. Let M^{y) denote the principal sub- 
matrix of M'^^{y) indexed by the subset {(Oa') | a' G N"}; in other words, 
one may think of M'^{y) as being indexed by N" with (a', /3')th entry ya'p'', 
let us call M'^{y) a pruned complex moment matrix. These three types of 
matrices M^{y) {K = R, C) and M^'^{y) will play a central role in our treat- 
ment. It will be convenient to think of M^[y) as being indexed by T„ and 
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of M^'^{y) as being indexed by 

fn := \a,a' e N"} C C[x,x]. 

Thus T„ ~ and we view x as a complex variable in the complex 

case. Recall that one says that '/ G K[x] lies in the kernel of M^(y)' 
if M^(y)Yec{f) = 0. Similarly, one may identify a polynomial {x,x) i— > 
f{x,x) = J2a,a' fa,a'X'^x°'' with its sequence of Coefficients vec(/) = {fa,a')a,a' 
which allows us to say that '/ € C[a;, x] lies in KerM2C(y)' if M2C(y)vec(/) = 
0. 

We also need truncated moment matrices. For an integer t > 0, M^{y) 
denotes the principal submatrix of M^{y) indexed by T„^t and Mf'^{y) de- 
notes the principal submatrix of M'^^{y) indexed by the set T„^( := {x°'x°' \ 
a, a' G N", |a| + \a'\ < ~ Tan,*- Given h G M[x], h{x) = iZf^hjixl^ , and 
y G M^", define hy G M^" by 

hy := M^{y)vec{h); that is, {hy)a = ^ hVcc+lS for a G N". 

Similarly, given h{x,x) = J2i3 /s' hpis'X^x^' G C[x,x] and y G C^^", define 
hy G C^'" by 

hy := M2^(y)vec(/i); that is, ihy)aa' = X hp'ya'+p,a+p' for a, a' e N". 

When /i G C[x] (i.e., /i/j/?/ = if /3 / 0), M^(y)vec{h) is the projection of 
M'^'^ {y)vec{h) onto the coordinates indexed by the pairs (a, a') with a = 0. 

3.2 Measures and kernels of moment matrices 

For a Hermitian matrix A, write A ^ if ^ is positive semidefinite, i.e., if 
u*Au > for all u G C" (or u G when A is real valued). 

The real case. For G C", set := {v")aeN" and := ('y")tteN" for 
an integer t > 0. Let ;U be a positive measure on M" with finite support; 
say, n = J2veW ^-"^^ where A„ > and W C M", \W\ < oo. The sequence 
of moments of the measure n is the sequence yf^ G M^" defined by {y'^)a '■= 
J x^dj^i = ^^gyj/ Xyv" for a G N"; (y'')o = YuveW ^^'^ total mass of the 
measure, equal to 1 if /x is a probability measure. We have 

vew 

Moreover, M^{y^^) = E^evK ^vCvCl t and 

KciM^{y^) = {/ G M[a;] | f{v) = e W} = I{W), 
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KerMf (y'') = I{W) D R[x]t 

(which follows from the fact that vec(/)^Mf (C2t,^,)vec(/) = /(v)^ for / G 
M[x]f). Given polynomials hi, . . . ,hm & let dj,d be defined as in (1.2) 
and, for t > d, set 

:= {y G M^^* \ yo = I, Mf{y) ^ 0, Mf_, (/i,y) = (j = 1, . . . , m)}. 

(3.2) 

Then, Kf- is a convex set which contains the vectors C2t,v for all v ^¥^{1). 
The following geometric observation, which indicates how the real radical 
ideal of / relates to the kernel of moment matrices, will play a central role 
in our approach. 

Lemma 3.1. Let I = {hi, ... , hm) ^ W\x\, t > d, and let y G Kf- for which 
rankMj^(y) is maximum. Then, KerM]^(y) C KeTM^{z) for all z G Kf. 
Moreover, KerMf (y) C I{Vth{I)). 

Proof. Let z G Kf. Then, y' := ^y + z) G Kf and KerMf (y') = 
KerMf (y)nKerMf (z) C KerMf (y). As rankMf (y) > rankMf (y'), equal- 
ity KerMf (y) n KerMf (z) = KerMf (y) holds, which implies KcrMf (y) C 
KeiMfiz). As C2t,v G Kf for ah v G V]r(/), this imphes Ker Mf(y) C 
n^gVi(/)KerM/^((^2t,D) which in turn is contained in I{Vk{I)). □ 

The complex case. Let ^ be a positive measure on with finite support; 
that is, ji = ^y^iY XySy where > and W C C", \W\ < oo. One can now 
define the doubly-indexed sequence of moments y'^ G C^^" of the measure jx 
by {y^)aa' '■= J x°'x'^ djjL = Y^y^w for ct, € N". Thus y^ satisfies 

(3.1) and 

Therefore, M^^(y^) = ^y^w '^vCv ® Cv {Cv Cv)'^ t and, in particular, 
M^iy^') = ^vCvC >1 0. Moreover, 

KerM^(y^) = {/ G C[x] \ f{v) = Vf G VF} = I{W), 

KevM'^^iy^') = {/ G C[x,x] \ f(v,v) = € W} 

(using the fact that vec(/)*M2C(^- (^)vcc(/) = |/(f,U)P for / G C[a;,^]). 
Given polynomials hi, ... , hm G C[x], t > d, define the sets 

:= {y G C^i" I yo = 1, (3.1), Mf(y) t 0, 

^f-d,(M = 0(j = l,...,m)}, 
Kf := {y G C^i" | yo = 1, (3.1), Mf{y) ^ 0, ^^"^^ 
M2_C^(M = 0(i = l,...,m)}. 

Hence, -ftT^"^^ CI K^' arc both convex sets. The following analogue of Lemma 3.1 
holds in the complex case; we omit the proof. 
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Lemma 3.2. Let I = {hi,... , hm) C C[x] and t > d. 

(i) Let y G for which TankMf{y) is maximum. Then, KerM^^(y) C 

KevMfiz) for all z e K^. Moreover, KevMfiy) C L{Vc{I)). 

(ii) Let y G Kf^ for which rankMf^(?/) IS maximum. Then, Kei M^^{y) C 

Ker Mi^{z) for all z G K^^. Moreover, Ker Mf{y) C I{Vc{I)). 

Link between the real and complex cases. As shown e.g. in [14] the 

complex moment problem in C" can be reduced to the real moment problem 
in M^**. Let us sketch the main idea. For a, a' G N", define the polynomials 

$("«')(.,^) := (^)"(^)" = i;4"A'-'^' 

f3(3' 

The following can be easily verified: For all x G C", and all (n, v) G M", 

In addition, $(""'+/^/^') = and = \I/(""')\i/(/3/3'). More- 

over, for every a, a' G N", 

Next, given y G C^^", define the linear mapping Ly : C[x,x] C by 
Lyif) = T,p/3' UP'VPP' foi" / ^ with f{x,x) = Y^pp, fpp'x'^xl^', and 

the mapping : " — ?/ 1— a := (^(y) with: 

a«„, = L,($("«')) = ^ {aa' G N^"). (3.4) 

Notice that a G M^^" if y satisfies (3.1). Conversely, given a G C^^", let 
La •.C[u,v\ — ^ C be the linear mapping 

9 (= '^9pp'U^v^') ^ La{9) ■= '^9pp'a(3p', 9 ^ C[u,v], 

and the linear mapping i/j : C'^^" —>■ C'^^", a t-^ y := ip{a) hy 

W = L„(^(""')) (aa' G N^"). 

Notice that y satisfies (3.1) whenever a is real valued. The mappings ip and 
i/j are inverse bijections between the set of sequences in C^^" satisfying (3.1) 
and R^^". Based on the above observations, we can now verify that 

M2C(y) >- M^{a) y 0. 
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Assume first M'^'^{y) y and let / G M^'" arbitrary. Then 

aa',l3/3' aa' 

= Ly{C£ = vec(5*)M2C(y)vec(5) > 0, 

aa' 

with g{x,x) := Eaa' /""'^^""'Ha^.^)- This shows that M^{a) h 0. Con- 
versely, assume M^{a) >: 0, and let g G C^^" arbitrary. Then 

g*M'^'^{y)g = ^ g^gp^'Va'+p^a+p' 

= Lai E W5/3/3'*^"'"+^^'^)=i^a(^^), 

with h{u,v) := ^^a' daoi'^^^^'Kujv). Now h = h\ + ih2 with /ii,/i2 € 
M[u, t;], and so 

L„(M) = La{hl + hl) 

= vec{hifM*{a)vec{hi) + vec(/i2)^M'*(a)vec(/i2) > 0, 

which shows that M'^'^{y) ^ 0. 

Finally, y is the sequence of moments of a measure on the set W C 
C" if and only if a is the sequence of moments of a measure on the set 
{(Im(i;),Re('y)) \ v e W} C M^n. (Use the fact that, if y = (g) (y, then 

« = C(Im(i;),Re(j;))-) 

3.3 Flat extensions and finite rank moment matrices 

Given a Hermitian matrix A and a principal submatrix B of A, one says that 
A is a flat extension of B if rankvl = rankZ?; then A y <^=^> 5^0. We 
begin with two fundamental results of Curto and Fialkow [12] about finite 
rank moment matrices, where this notion of flat extension plays a central 
role. See [27] for a short proof of Theorem 3.3 and [29] for an exposition of 
Theorem 3.4. 

Theorem 3.3. (i) IfM^{y) >z and rankM'*(y) < oo, theny = Ei,eiy 

for some finite set C M" and Xy > 0, \W\ = rankM'^(y), and 
KerM^iy) = I{W). 

(ii) If M^^iy) h and TankM^^{y) < oo, then y = Ev&w Cv ^ Cv 
for some finite set W Q and Xy > 0, \W\ = rankM^^(?/), and 
KevM^iy) = I{W). 

Theorem 3.4. (i) If Mf{y) y and rankMf(y) = TankM^_^{y), then y 
can be extended in a unique way toye M^" such that M^{y) is a flat 
extension of Mf{y) (and thus M^{y) yO). 
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(ii) If Mf'^{y) >i and rankMj^'^(y) = rankM(^5i(y), then y can he extended 
m a unique way toy e such that M^^{y) is a flat extension of 

Mf^{y) (and thus M^^{y) ^0). 

The following lemma taken from [12] shows that the kernel of a truncated 
moment matrix enjoys ideal-like properties. 

Lemma 3.5. (i) Let Mf{y) yO,f,ge R[x], with deg{fg) <t-l. Then, 
(y)vec(/) = ^ Mf{y)Yec{fg) = 0. 

(ii) Let Ml^{y) yO,f,ge C[x,x], with deg{fg) <t-l. Then, 
M^^{y)Yec{f) = ^ Mf{y)Yec{fg) = 0. 

Proof. Set h := fg. (i) As deg(/i) < t - 1 and Mf{y) y 0, it suffices to 

show Mil {y)vec{h) = 0. Moreover it suffices to show the result for g = xf, 
in this latter case one can verify that, for a G ^^-n (-^t^i(j/)vec(/i))Q = 
(Mf (y)vec(/))«+e, = 0- _ 

(ii) Similarly assume g = Xi or Xi. For aa' G Tn^t-i, {M^Si{y)^Qc{h))aa' 
is equal to {M^^{y)vec{f))a+e, a' if 5 = Xi and to {M^^{y)Yec{f))a a'+a if 
g = Xi, thus to in both cases. □ 

Proposition 3.6. KcrM'*(y) is an ideal in M.[x], which is real radical if 
M^{y) y 0. Assume M^{y) t and rankM'^(y) = TSinkM^_^{y) for some 
integer t>l. Then, KerM'^(y) = (KerMf (y)) and, for B QT^, B indexes 
a (maximum) nonsingular principal submatrix of M*{y) if and only if B is 
a (maximum) linearly independent subset oflS.[x]/KevM*{y). 

Proof. We use the (easy to verify) identity: 

Yec{h)'^ M^{y)vccipq) = vec{hq)^ M^{y)vec{p) 

for p,q,h £ R[x]. Up £ KerM^{y), q £ R[x], then vec(/i)^M'^(y)vcc(pg') = 
vec(/ig)"^M'^(j/)vec(p) = for all h £ R[x], which implies M^{y)Yec{pq) = 
and thus pq £ Ker M^{y). This shows that KerM'^(y) is an ideal. As- 
sume now M*{y) >: 0; we show that KerM'^(y) is real radical. In view 
of Lemma 2.1, it suffices to show that if Yli=iPi ^ KcrM'^(y) for some 
Pi £ R[x], then pi £ KeiM^{y). Indeed, = vec(l)^iV#(y)vec(X;iLi Pi ) = 
^^^j^ vec(pi)-^M''^(y)vec(pi) implies vec(pi)'^M'^(y)vec(j3i) = and thus £ 
Kei M^{y) for aU i. 

Assume rankM^{y) = rankMf^_^{y) =: r and set J := (KerMf(?/)). 
Obviously, J C KerM''^(y); we show equality. For this, let B C Tn,t-i index 
an r X r nonsingular principal submatrix of M^{y). We show that, for all 
a G N", G Span]g(B) -|- J, using induction on |a|. This holds for |a| < t 
by the definition of B. Assume |a| > t + 1 and write = XiX^ . By the 
induction assumption, x^ = J2xf^eB^f^^^ ^ where q £ J, Cfs £ R. Thus, 

= X]a;/3eB C/jXix'^ + Xiq. Here, Xiq £ J and XiX^ G Spani[j(i3) -|- J since 
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deg{xixf^) < t, which impUcs G Span]g(;B) + J. Thus we have shown that 
M[x] = SpanK(.S) + J. As KerM'*(y) n Span]R(;B) = {0}, this impUes easily 
that KerM^(y) = J. 

For H C T„, it is obvious that B indexes a nonsingular submatrix of 
M^(y) if and only if B is linearly independent in ]R[x]/KerM'^(|/). The 
last statement of the lemma now follows since dim]R[x]/KerM'^(y) = r (as 
KerM'*(y) is radical and using the identity \Vc{KQTM*-{y))\ = rankM'^(y) 
from Theorem 3.3). □ 

Proposition 3.7. KeTM'^'^{y) is an ideal in C[x,x]. If M'^^{y) y 0, then 
KerM2C(y) 

is a radical ideal in C[x, x] and thus KerM'^(y) is a radical ideal 
in C[x]. Assume, moreover, rankM'^(y) = rankMj^^(y) for some integer 
t > I. Then, KerM'^(y) = (KerMf (y)) and, for B C T^, B indexes a 
(maximum) nonsingular principal submatrix of M^{y) if and only if B is a 
(maximum) linearly independent subset of C[x]/KeTM'^{y). 

Proof. As in the real case, we use the following (easy to verify) identities: 
For h,p,q e C[x,x], 

vec(/i)*M2C(y)vec(pg) = Yec{hp)* M^^ {y)Yec{q) , 
vec(p)*M2C(y)vec(^?) = vec(l)*M2C(y)vec(^?^), 
vec(p^)*M^'^(y)vec(p^) = vec(pp)*M^'^(y)vec(pp), 

where p E C[x,x] is defined as p{x,x) := p{x,x). This implies directly that 
KerM^'^(j/) is an ideal. Assume now M'^'^{y) ^ 0; we show that KerM^'^(y) 
is radical. In view of Lemma 2.1, this follows from the following fact: 

p2 G KerM2C(y) 

= vec(p2)*Af2C(y)vec(p2) = Yec{pp)* M'^^ {y)\ec{pp) 
pp G KerM2^(7/) 

= vec(l)*M2C(y)vec(j9^) = Yec{p)* M"^^ {y)Yec{p) 
^pGKerM2C(y). 

The proof for the last statements of the proposition is identical to the proof 
of the corresponding statements in Proposition 3.6 for the real case. □ 

Without the assumption M^'^{y) y 0, KerM'^(y) is not necessarily an 
ideal in C[x]. Indeed, for the sequence y G " defined by yaa' '■= 
if a = or a' = 0, and y^a' ■= 1 otherwise, M^{y) >z 0, M'^^{y) ^ 
and KerM*'(y) is not an ideal (e.g., 1 G KerM'^(y) while any nonconstant 
monomial does not lie in KerM*'-'(y)). We mention for further reference the 
following corollary, and wc conclude the section with the proof of Proposition 
1.1 and a lemma about the '(real) radical'-like property of the kernel of a 
positive semidefinite truncated moment matrix. 

Corollary 3.8. (i) Assum,e Mf{y) h andrankMf{y) = rankMf_^{y) =: 
r. Then, J := (KerM^(y)) is real radical and zero- dimensional, 
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dimM[a;]/J = r, J n R[x]s = KerMf (?/) and, for B C T„,^, B in- 
dexes a (maximum) nonsingular principal submatrix of Mf-{y) ■^=^ B 
is (maximum) linear independent in M.[x]/J. 

(ii) Assum,eM^^{y) >z andiankM^^{y) = rankMf_i(y) =: r. Then, J : = 
(KerMf(y)) is radical, dimC[x]/J = r, Jr\C[x]s = KerMf(?/) and, 
forB C Tn,s, B indexes a (maximum) nonsingular principal submatrix 
of M^{y) <^=^ B is (maximum) linearly independent inC[x]/J. 

Proof. We prove only (i). By Theorem 3.4, y has an extension y G M^" 
such that M^{y) is a flat extension of M^{y). By Proposition 3.6, the 
ideal KerAf^(y) = (KerMs(y)) =: J is real radical and zero-dimensional, 
dimM[x]/J = r, and J n R[x]s = KeiM^iy) D R[x]s = KerMf (y) = 
KerMf(y). □ 

Proof of Proposition 1.1. By Proposition 3.6, J := KerM'^(y) is a real 
radical ideal, since M^{y) y 0. As = M^{hjy) = M^{y)vec{hj) for all 
j, we have I ^ J, which implies that Vir(J) C is finite. As J is real 

radical, we deduce that Vc(J) = V^{ J) C W^. Hence J is zero-dimensional 
and /(Vk(/)) C J since Vc{J) C V^{I). Set r := dimM[x]/J = |Vc(J)| < 
|Vm(I)|. Let i3 C T„ be a linear basis of M[x]/J, \B\ = r. Then the columns 
of M^{y) indexed by B form a basis of the column space of M^{y) and thus 
rankM''^(y) = r. Moreover, r = |11r(/)| if and only if Vc(J) = Vr{I) which 
in turn is equivalent to J = I{Vr{I)). Now, this maximum rank |Vr(/)| is 
reached by the sequence y := y^ = ^y^Vu{i) '^ith > which indeed 
satisfies (1.1). □ 

Finally, it is useful to observe that the kernel of a positive semideflnite 
truncated moment matrix enjoys the following '(real) radical'-like property. 
We omit the proof whose details are straightforward. 

Lemma 3.9. (i) Assume Mf (y) >z and let p, qj G R[x], f := Qj 
with m G N, m > 1. Then, f G KerMf (y) =^ p G KerMf (y). 

(ii) Assume M^'^{y) >z and let p G C[x], m e N, m > 1. Then, p"^ G 
KerM2C(y) ^p^ KeiM'^^{y). 

4 A semideflnite characterization of the (real) rad- 
ical ideal via moment matrices 

In this section we present a semideflnite characterization of the real radical 
ideal /(Vk(/)) of an ideal / C M.[x\, as well as a numerical algorithm for 
computing a set of generators. It turns out that the method also applies 
to the radical ideal /(Vc(/)). Our strategy is to obtain /(Vic:(/)) (K = R 
or C) as the ideal generated by the kernel of some suitable moment matrix 
M^{y) where y G K^. Sections 4.1-4.3 contain some results ensuring that 
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the moment matrix M^{y) has the desirable properties for achieving this 
task and Section 4.4 describes our algorithm. 

4.1 Weakest set of conditions 

Throughout, / = {hi, . . . , hm) is an ideal in for which we want to find 
the radical ideal I{V^{I)), K = R or C. Recall the definition of d in (1.2). 

Proposition 4.1. Let t>d, l<s<t, for which rankM^(y) is 

maximum, and B C Tn,s-i index a maximum nonsingular principal subma- 
trix of Mf_-^(y) with border dB defined as in (2.3). Assume (i)-(iii) below 
hold: 

(i) B is an order ideal. 

(ii) The principal submatrix of M^{y) indexed by BUdB has the same rank 

as M^_-^{y); that is, with B := {bi, . . . , 6;v} and dB := {ci, . . . , ch}, 
there exists a polynomial gj G Ker M^{y) of the form gj{x) = cj — 
^iLi aijbi (i.e., G := {gi, . . . ,gH} is a B-border prebasis). 

(iii) The formal multiplication matrices Xi, . . . ,Xn defined from G commute 
pairwise. 

Then G is a border basis of J := {G) C /(Vk(/)), B is a linear basis of 
M.[x]/J, and one can extract (using the formal multiplication matrices) the 
set W := VciJ) which satisfies Vk{I) C W and \W\ < iankMf_^(y). More- 
over, if \Vk{I)\ = \W\ = rankMf_i(y), then Vk{I) = W, I{Vk{I)) = J, and 
G is a B-border basis of I{Vk{I)). 

Proof. Theorem 2.8 gives directly that G is a border basis of the ideal 

J = {G) and that ;B is a linear basis of M[x]/J. Moreover, the matrices 
Xi, . . . , Xn coincide with the multiplication matrices in ]R[x]/J w.r.t. the 
basis B. Thus one can compute the set W = Vc(J) from their eigenvectors. 
By construction, J C (KerMf (y)) C (KerMf (y)) C I{Vk{I)), where the 
last inclusion follows from Lemmas 3.1 and 3.2 (i). This implies Vk{I) Q 
Vc{J) = W. Moreover, \W\ < dimK[x]/J = \B\ = rankUft^iy). 

U \ Vk{I)\ = \W\ = rankM|ii(y), then W = Vk{I) and J is radical since 
dimK[a;]/J = \Vc{J)\, which implies I{Vk{I)) = I(W) = I{Vc{J)) = J. D 

In the next two subsections we give simple rank conditions (4.1), (4.5), 
which ensure that the conditions of Proposition 4.1 hold. (See Remark 4.13 
for details.) For the sake of clarity, we treat the real and complex cases 
separately. 

4.2 Chctracterizing and computing the real radical /(Vk(/)) 
Assume I = {hi, . . . , km) is an ideal in M[x], with hi, ... , km G M[a;]. 
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Proposition 4.2. Let t>d and y G Kf- for which rankMj^(y) i 
If, for some 1 < s <t, 



IS maximum. 



rankMf (y) = rankMf_i(y) =: r 



(4.1) 



then /(V]r(/)) D (KerMf (y)) =: J. One can compute the set W := Vc{J) 
which satisfies C W and \W\ = r. Moreover, if Vr{I) = W then 

I{Vr{I)) = J. 

Proof. By Lemma 3.1, KerMf(y) C KerMf (y) C implying 
J Q and thus Vr{I) C W. By Corollary 3.8 (i), as J is radical, 

\W\ = dimR[x]/J = r, and Vr{I) = W implies I{Vr{I)) = I{W) = J. □ 

One can verify (using Lemma 3.5) that, if (4.1) holds for some s < t — 2, 
then it also holds for s = t — 1. Hence it suffices to check whether (4.1) 
holds ioT s = t — 1 or t. In Lemma 4.3 below, we observe that, if assumption 
(ii) in Proposition 4.1 holds for s < t — 1, then in fact (4.1) holds and thus 
Proposition 4.2 applies. However, it may be that Proposition 4.1 applies to 
the case s = t while (4.1) does not hold; see Example 5.3 (for relaxation 
order t = 2) for such an instance. In Remark 4.13 below we see that the 
converse of the next lemma holds. 

Lemma 4.3. In Proposition 4-1, if assumption (ii) holds for s <t — l then 
TankMf_^{y) = rankMf (y). 

Proof. We have to show that Tn,s Q SpsLii^^B) + KerMs(y). By the 
definition of B, any G Tn,s-i lies in Span]g(H) + KerM5_i(y). For G 
Tn,s, write a;° = xix^ where x^ G Tn,s-i- Thus x^ = Xlz/^eB '^/^^'^ ~^ P 
where p G Kei Mf^_-^^{y), C/j G K. Therefore, = J2xi^<^B ^1^^'^^^ + xip. As 
xix^ G BVJdB, assumption (ii) implies that xix^^ G Spanjj(;S) + KerM|^(y). 
As deg(2:ip) < s < t — 1, it follows from Lemma 3.5 that xip G KerM|^(y). 
Therefore, x'* G SpanK(H) + KerMs(y). □ 

By strengthening the rank condition (4.1), one can show that J = 
/(Vm(/)), i.e., the desired real radical ideal is found. 

Proposition 4.4. Let t>d and y G Kf for which rankMj*(y) is maximum. 
Assume that, either (4-1) holds for some 2d < s < t, or 



for some d< s <t. Then I{Vr{I)) = (KerMf (y)) (and one can find Fr(/) ). 
Moreover, rankMf (y) = \Vr{I)\. 

Proof. In view of Proposition 4.2, there remains only to show the in- 
clusion I{Vr{I)) ^ J := (KerM^(y)) or, equivalently (since J is radical), 
W ■.= Vc{J) CVu{I). We already know that C M" since J is real radical 



rankMf (y) = rankMf_rf(y) 



(4.2) 
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and zero-dimensional (see Corollary 3.8). We now show that W C Vc{I)- 
Assume first that (4.1) holds for s > 2d. As Mf'_j,\hjy) = 0, we have 
{hjy)a = for all |a| < 2t — 2dj and thus for all |a| < 2dj. Hence, 
(y)'^ec(/ij) = and thus hj E KerMf(y). Therefore, / C J, giving 
W C Vc{I). 

Assume now that (4.2) holds for sonic d < s < t. Let (v G W) be 
interpolation polynomials, i.e., Pv{w) = dv,w for v,w £ W. As observed in 
[29, Lemma 27], one can assume that deg(p^) < s—d. (Indeed, let B C T^^s-d 
index a maximum nonsingular submatrix of Mf-{y); then B is a basis of 
M[x]/ J by Corollary 3.8 and one can replace py by its residue modulo J w.r.t. 
B.) ^From Theorems 3.3, 3.4, we know that iya)aeT„,2s = Y^veW ^vC2s,v 
where > 0. Hence, = vec{py)^ Mf_^ {hjy)vec{py) = hj{v)Xy implies 
hj{v) = for all j and thus v G Vc{I), which shows W CVc{I)- □ 

We now formulate an analogous result for the ideal /(Vk(/) n S), where 
5 := {x G M" I hm+i{x) > 0, . . . , hm+k{x) > 0} is a semialgebraic set, with 
hm+i, ■■■ , hjn+k S ^[x]. For this define the set 

K^S ■■= Kf n {y I Mt_d^ {hjy) ^ {j = m + 1, . . . ,m + k)} (4.3) 

foi t > d := max dj. 

j=l,...,m+k 

Proposition 4.5. Let t > d and y G Kf'g for which TankM^{y) is maxi- 
mum. 

(i) Assume (4.I) holds for some 1 < s < t. Then J := (KerMf(j/)) C 

I{Vr{I) n S) and W := Vc{J) ^ Vu{I) n S with \W\ = rankMf (y); 
moreover, J = I{Vu{I) nS) ifW = Vr{I) n S. 

(ii) Assume (4.2) holds for somed<s<t. Then I{V^{I)r\S) = {KevMf{y)) . 

Proof. (i) The inclusion KerMf (y) C /(Vs;(/) n S) follows from the 
maximality of the rank of Mf'{y) and the fact that C,2t,v G ^f's ^ 
V G Vk{I) n S. This gives J C I{Vk{I) D S) and thus W D Vk{I) D S. 
Equality W = Vr{I) n S implies I{Vr{I) n S) = I{W) = J (as J radical). 
This concludes the proof of (i) . The proof for (ii) is analogous to that of the 
corresponding statement in Proposition 4.4. □ 

To conclude we show that, when Vm(7) is finite, then condition (4.2) 

is satisfied for t large enough. That is, the conclusion of Propositions 4.4, 
4.5 holds: the real radical ideal /(Fm(/)) or I{Vr{I) n S) = (KerMf (y)) is 
found. 

Proposition 4.6. Assume |Vr(/)| < 00. 

(i) //Mr(/) = then K^g = ^ fort large enough. 
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(ii) If Vm.{I) 7^ then, for t large enough, there exists d < s <t such that 
rankMf (y) = ia.nkMf_^{y) for all y e K^g. 

Proof. Assume t >2d and let y G Kfs- Then, as observed in the proof 

of Proposition 4.4, hi, ... , hm G KerM]^(y). We first show that, for t large 
enough, KerM]*(y) also contains a given basis of the ideal /(Vk(/)). 

Claim 4.7. Let {gi, ■ ■ ■ , gk} be a basis of the ideal /(Vk(/)). There exists 
to G N such that gi,. . . ,gk G KcrMf (y) for all t>tQ. 

Proof. Let / G {1, . . . , fc}. By the Real Nullstellensatz, there exist nii G N, 
mi>\ and polynomials ai, uj^ {j < m) for which gf"^^ + ai = Xljli ^j'^^i 
and (7/ is a sum of squares. Set to := 1 + maxi<fcj<m(2d, deg(g';^™'), deg(a";), 
deg{uf^hj)) and let t > to. As deg{uf''hj) < t - 1 and hj G KerMf (y), then 

ufhj G KerMf(y) by Lemma 3.5. Hence, gf""^ + ai e KerMf(y) which, 
using Lemma 3.9 (i), implies gi G KeTMf'{y). □ 

If Vr{I) = 0, then {1} is a basis of /(FrW) = ^N- Hence 1 G 

KerM]'^(y), implying yg = 0, which contradicts the fact that yo = 1 for 
y G -fCf^. Therefore, i^f^, = for t > to, which shows (i). 

Assume now Vm.{I) ^ 0. Let {gi, . . . , gk} be a Grobncr basis of I{Vr{I)) 
for a graded monomial ordering ordering and let B be the correspond- 
ing set of standard monomials. Thus B is a basis of W[x]/I{Vr{I)); set 
de := maxftggdeg(6) (which is well defined as S / 0). We can write any 
monomial as x" = r^") + ^1^1=1 Pi"'^ 9h where r^") G Span]g(6), G W[x] 
and dcg{p^"^ gi) < dcg(x"). Set ti := max((i0 + ti, to) and let t > ti + 1. 
Consider a G T„,(,. As dcg{p^i^^ gi) < ti < t- 1 and gi e KerMf (y), we have 

pl^'^gi G KerMt(y) and thus - r^") G KerMt(y). As deg(r(")) < de < 
ti — d, this shows that the ath column of M^{y) is a linear combination 
of columns indexed by Tn,ti- Therefore, rankM('^(y) = rankM]^_^(y), thus 
proving (ii). □ 

Remark 4.8. Detecting existence of real solutions. Hence one can 
detect the existence of real solutions via the following criterion: 

Vu{I) = <s=^ iff = for some t. (4.4) 

(The 'only if part follows directly from Proposition 4.6 (i), while the 'if part' 
follows from the fact that C2t,v G Kf' for any v G Vr{I).) Note moreover 
that, when Vk{I) = 0, none of the flat conditions (4.2), or (4.1) with s > 
2d, can hold; indeed, under either of these two conditions, one would have 
|Vm(/)| = rankM^(y) > 1 by Proposition 4.4. Consider as an illustration 
the following small example: I = {h := xf + 1) C ]R[xi] with Vr{I) = 0. 
For t > 1, if y G i^t then Afi(y) >z implies y2ei > 0, while (/iy)o = gives 
2/2ei + 1 = 0, yielding a contradiction. Hence, Kf = ^ for any t > 1. 
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Remark 4.9. Proposition 4.6 remains valid under the weaker assumption 
n S"! < cx) if, in the definition of the set Kf-g in (4.3), we add the 

constraints Mt-dSVeV) ^ for e G {0, l}'^, after setting := Y[i=i Kn+i- 
The proof is analogous, except we now prove in Claim 4.7 that KerMj^(j/) 
contains a given basis of the ideal 7(Vk(7) 5). To show this, instead of 
the Real Nullstellensatz, we now use the Positivstellensatz (see Stengle [39] ) 
which in our case can be formulated in the following way: For g G M.[x\, 
g G /(Mr(/) n S) if and only if -g'^'' = Yj'jLi Ujhj + I]ee{o,i}fe ^ePe for some 
r G N \ {0}, Uj,ae G R[x], with ag s.o.s. 

4.3 Chctracterizing and computing the radical I{Vc{I)) 

Using complex moment matrices, we can formulate analogues of Proposi- 
tions 4.2, 4.4, 4.6 for the radical ideal I{Vc{I)); the proofs of the first two 
results being similar are omitted. 

Proposition 4.10. Let t > d and y G Kf^ for which TaxikM^'^{y) is maxi- 
mum. If, for some integer 1 < s < t, 

rankMf'^(y) = rankMf_i(y) =: r (4.5) 

then I{Vc{I)) ^ (KerMf (y)) =: J. One can compute the set W := Vc{J) 2 
Vc{I) which satisfies Vc{I) C W and \W\ = r. Moreover, if W = Vc{I) 
then I{Vc{I)) = J. 

Proposition 4.11. Let t > d and y G Kf^ for which rankMj^'^(y) is maxi- 
mum. Assume that, either (4-5) holds for some 2d < s < t, or 

rankMf^(y) = rankMf_^(y) (4.6) 

for some d < s < t. Then I{Vc{I)) = (KerM^(y)) (and one can find Vc{I)). 

Proposition 4.12. Assume |Vc(/)| < oo. 

(i) // Vc{I) = 0; then Kf^ = ^ for t large enough. 

(ii) If Vc{I) 9 then, for t large enough, there exists d < s <t such that 

rankM|C(y) = rankMf_^(y) for all y G Kf^. 

Proof. Let t > 2d and y G K^^. Then, hi,..., km G KerMf (y), since 

{Mf{y)hj)a' = {hjy)oa' = for \a'\ < t < 2t — 2dj by the assumption 
^i-dji^jV) = 0. Hence, hj G KerM^^C^^) ^j^^g 75- ^ KerM^'^{y) too. 

Let {gi, . . . ,gk} be a Grobner basis of I(Vc(-^)) for a graded monomial 
ordering. Analogously to Claim 4.7, one can show, using Hilbert's Null- 
stellensatz, the existence of G for which gi G KerM^^'^(y) for all I 
and t > to. If Veil) = 0, then 1 G KcrM^^{y) which implies yo = 0, 
thus showing K^'^ = 0. Assume now that Vc(/) / 0. Let B be the 
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basis of C[x]/I(Vc{I)) for the ehosen monomial ordering and set djs := 
maxfeeB deg(6) (which is well defined as ;S / 0) and C := {W \ b,b' e B}. 
Then any monomial x".x"' can be written = r^""') + Ylu'=i ^L"" ^W9i' 

where ■* G C[x,x], r^"" ^ G Spanc(C), and deg(tt|;?" ^ViQi') ^ l"^! + l^^'l- 
Let ti := max{3d,d + 2dB) and t>ti + l. Then x V' - r(""') G KevM^'^iy) 
whenever |q; + a'| < ti which, together with deg(r'^""^ < 2dB < h — d, 
shows that rankMj^^'^(y) = rankMj^^'^^(y) =: r. 

There remains to show that rankM^^_^(y) = r. Applying Theorems 3.3, 
3.4, there exists W C C^, \W\ = r, Xy > (v £ W) such that, if we set 
y ■= Y^vew^vCv^ Cv, then M'^'^{y) is a flat extension of M^^{y). This 
implies M^^_^{y) = J2vsW ^vCti-d,vCl-d,v G KcrM^^_^(y) (since 

ti — d > deg{hj) as ti > 3d), we deduce that hj(v) = for all f G 
and thus W C Vc(/). We now show that the vectors Cti-d,v [v G W) are 
linearly independent, which implies that rankMj^_^(y) = \W\ = r, thus 
concluding the proof. For this, consider interpolation polynomials G C[x] 
{v G W), i.e., satisfying Pv{w) = 6^^^ for v,w € W. One may assume that 
deg(pt;) < ds (replacing if necessary by its residue modulo /(Vc(-^)) with 
respect to the basis B). Assume J2veW '^■"iti-d,v = for some G C; we 
show that all c^'s are zero. As ti — d > c/g, we can take the scalar product 
with vec{pu,) {w G W) which yields = I^^g^F (^Pw{v) = Cw □ 

Remark 4.13. Under condition (4.1) or (4.5), the assumptions (i)-(iii) of 
Proposition 4.1 hold. Namely, one can construct an order ideal B C T„ s_i 
indexing a maximum nonsingular principal submatrix of AI^_i{y). More- 
over, one can choose B = By, the set of standard monomials for the ideal 
J := (KerM|^(y)) with respect to a graded lexicographic order; this is pos- 
sible since J is zero-dimensional and there is a basis in T^^^-i for ]fC[x]/J, 
which implies By C Tn,s~i by Lemma 2.5. Such a basis By can be found 
using the greedy sieve algorithm described in Section 2.4. The execution 
of the algorithm requires checking whether some set T C Tn,s-i is linearly 
independent in ]K[a;]/J. Using Corollary 3.8, this can be checked by testing 
whether T indexes a nonsingular principal submatrix of Mf'{y), thus by a 
rank computation on Mf'iy). Finally the formal multiplication matrices as 
defined in Proposition 4.1 coincide with the (usual) multiplication matrices 
in ]K[x]/J and thus they commute pairwise. 

Therefore, the conclusion of Proposition 4.1 also applies under (4.1) or 
(4.5): liW = Vk{I) then one can construct a border basis of J = /(Vk(-f))- 
Moreover, when using By, one can also construct the reduced Grobner basis 
of J for the graded lexicographic monomial ordering. A sufficient condition 
for W = V-k{I) is given above in Proposition 4.4 (K = M) and Proposi- 
tion 4.11 (K = C). 

In fact, as will be explained in Section 4.4.5, we can adapt this strategy 
to find a Grobner basis for an arbitrary monomial ordering. 
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Remark 4.14. The above results involve the matrix Ml'^{y) where the 
argument y G C^^"' is a complex sequence satisfying (3.1). As explained in 
Section 3.2 above, one may work instead with the moment matrix M^{a) 
where a G RN'" is the real sequence defined as in (3.4), and Propositions 
4.10 and 4.11 could be reformulated in terms of real sequences only. 

In fact, when the ideal / is generated by real polynomials hi, . . . ,hm, 
its set Vc{I) of complex roots is closed under complex conjugations, i.e., 
V e Veil) <S=^ V e Veil) and, for a polynomial / G C[a;], / G /(Vc(^)) if and 
only if its real and imaginary parts belong to /(Vc(/)); that is, it suffices to 
determine I(Vc(/))n]R[x]. For this, it suffices to consider real valued matrices 
Mfiy) (or Ml^iy)), i.e., with y G nM^'" (or K^CnM^ir) in Propositions 
4.1, 4.10 and 4.11. Indeed, one may e.g. easily verify that Lemma 3.2 
remains valid within the context of real polynomials and replacing Kf' (or 
K^^) by n M^2" (or K'i^ n M^^r). (Use here the fact that, since Veil) 
is closed under conjugation, then \iC,2t,v ® C,2t,v + C,2t,v ® (,2t,v) belongs to 

i^^cnRN^r.) 

4.4 Algorithm and implementation 

With the results of Sections 4.1-4.3 we have all the ingredients needed to 
compute the radical ideals /(Vk(/)) and /(Vc(-f)) of an ideal I given by its 
generators. We now describe the algorithm in more detail. 

For convenience, let Kt (resp., Mtiy)) stand for K^, Kf, Kf^ (resp., 
Mfiy), Mf{y), M^^iy)). For the task of computing liVuil)), we will use 
Kt = Kf' (and apply Propositions 4.1, 4.2, 4.4) and for the task of computing 
/(Vc(/)) we use Kt = Kf (and apply Proposition 4.1) or Kt = Kf^ (and 
apply Propositions 4.10, 4.11). The algorithm consists of five main parts: 
For a given order t> d, 

(i) Find an element y G Kt maximizing the rank of Mtiy). 

(ii) Check the ranks of the principal submatrices of Mtiy). 

(iii) Compute a basis for the column space of Mg^iiy) and the quotient 
space K[a;]/J (J = (KerMs(y)), for suitable I < s <t). 

(iv) Compute the formal multiplication matrices. 

(v) Construct a basis for the ideal J. 

In step (ii) we search for a submatrix Mgiy) of Mtiy) satisfying Proposi- 
tion 4.1 (i)-(iii), or the rank condition (4.1) (resp. (4.5)), or (4.2) (resp. 
(4.6)). Depending on what condition is satisfied, the algorithm returns a 
subideal J C /(Vk(/)) together with a superset W 5 V^il), or the desired 
radical ideal /(Vk(-^)) and the desired set of roots Mk(/). One can any- 
way verify a posteriori whether W = !%(/), simply by checking whether 
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hj{v) = for all j < m and v e W. In the sequel of this section we give 
more details about these different tasks. 

4.4.1 Finding y G Kt maximizing the rank of Mt{y) 

This first task can be cast as the problem of finding a feasible solution of a 
semidefinite program, that has maximum rank. For details on the theory and 
applications of semidefinite programming the interested reader is referred, 
e.g., to [43], [46]. It is a known geometric property of semidefinite programs 
that a feasible solution has maximum rank if and only if it lies in the relative 
interior of the feasible region and that such point can be found with interior- 
point algorithms using self-dual embedding (see, e.g., [15], [46]). Let us give 
some details. 

Consider a general instance of semidefinite program 

m m 

p* := inf ^ bjVj s.t. ^ yjAj -0^0 (4.7) 

and its dual semidefinite program 

d* := supTr(CX) s.t. Tr{AjX) = bj {j = 1, . . . ,m), X y 0. (4.8) 

Here, Aj,C,X are Hermitian matrices, b,y e M"*, X, y are the variables. 
Obviously, d* < p* (weak duality). There is no duality gap (i.e., p* = d*), 
e.g., when (4.7) is strictly feasible (i.e., 3y £ R™ with YJjLiVjAj - C >- ^) 
or when (4.8) is strictly feasible (i.e., 3X y feasible for (4.8)). When (4.8) 
is strictly feasible and d* < oo, then (4.7) attains its minimum, i.e., the set 
of optimal solutions is nonempty. The feasible region to (4.7) is the convex 
set 

K = {y\ Y:T=i Vj^j - C ^ 0} = {y I n*(Ef=i VjAj - C)u > Vn e K^}. 

Therefore, for y G K, y lies in the relative interior of K if and only if 
Ker(X;^i VjAj - C) C Ker(X;!^i zjAj - C) for all z G K or, equiva- 
lently, if Yl%^i Uj^j ~ CJ has maximum possible rank (same argument as 
for Lemma 3.1). 

Semidefinite programs can be solved in polynomial time to an arbitrary 
precision using, e.g., the ellipsoid method, whose running time is however 
prohibitively high in practice. Interior-point methods are now the method of 
choice for solving semidefinite programs. Assuming strict feasibility of (4.7) 
and (4.8), interior-point algorithms construct sequences of points on the 
so-called central path, which has the property of converging to an optimum 
solution of maximum rank [20]. One can also find a maximum rank optimum 
solution under the weaker assumption that (4.7), (4.8) are feasible (but not 
necessarily strictly feasible), if p* is attained, and p* = d* < oo. Indeed 
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one can then construct the so-called extended self-dual embedding which is 

a strictly feasible semidefinite program with the property that a maximum 
rank optimum solution to it yields a maximum rank optimum solution to 
the original problem (4.7) (see e.g. [15, Ch. 4], [46, Ch. 5]). 

For our problem of finding y G Kt maximizing rankMi(y), consider the 
semidefinite program 

p* := min 1 s.t. Mt{y) h 0, Mt_d. (h^y) = (j = 1, . . . , m), yo = 1, (4.9) 

where we add the condition (3.1) in the complex case. One can interpret 
(see, e.g., [25]) the dual of (4.9) as 

d* := max A s.t. 1 — A = s + YlJLi '^'j^'J ^'^^^ 1j polynomials 

deg(s), deg(gj/ij) < 2t, s is s.o.s. 

(4.10) 

where 's is s.o.s.' means that s can be written as a sum of squares, i.e., 
s = J2h l^/iP some polynomials Uh € K[x] or C[x,x]. Obviously, (4.9) is 
feasible if Mk(-^) / 0- Moreover, (4.10) is feasible (e.g. with A = 1, s = = 
as feasible solution) and, if Kf ^ 0, then p* = 1 is attained by the whole 
set Kf and p* = d* = 1. Hence an interior-point algorithm implementing 
the self-dual embedding technique applied to problem (4.9) is guaranteed to 
return the following information^: Either (i) y £ Kf maximizing rankMt(y), 
or (ii) a certificate that (4.9) is infeasible thus implying Vk{I) = 0- For our 
computations we use the semidefinite programming solver SeDuMi-1.05 [41, 
42] which has this feature. Practically, this means that the solution returned 
by the algorithm is very close to a maximum rank optimum solution. 
Remark 4.15. When using a semidefinite programming solver without 
the maximum rank property, one can recover a maximum rank solution 
to (4.9) from a feasible solution y to (4.9), using the following simple iter- 
ative algorithm. Let ui,. . . ,Up be a set of vectors that span KcrMt(y), set 
1 UiU*, and consider the semidefinite program: max(C, Mt{y)) sub- 
ject to y satisfying the constraints of (4.9). If the optimum value is equal to 
0, then y is in fact a solution of maximum rank. Otherwise, let yi be the op- 
timum solution returned by the solver; then KerMf(y) ^ KerMt{yi). Then, 
2/2 := ^{y+yi) is feasible for (4.9) andKerMt(y2) = KcrM^ (y)nKerMt(yi) C 
KerM((y). Hence we have found a feasible solution y2 to (4.9) for which the 
rank of Mt{y2) is larger than that Mt{y). Iterate replacing y by y2- 

4.4.2 Checking ranks of submatrices of Mt{y) 

Once a maximum rank matrix Mt{y) is found, one has to check if for some 
1 < s < t the conditions of Proposition 4.1 (i)-(iii) hold, or if (4.1) (resp. 

^Three options (I), (II), (III) are described in [46, Ch. 5, p. 119]; (i) corresponds to (I) 
and (ii) to (II), (III). Indeed, under (III) a certificate is reported that no complementary 

pair exists which imphes, in our case, that (4.9) is infeasible, since any y £ Kt together 
with the solution A = 1, s = Qj = to (4.10) makes a complementary pair. 
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(4.5)) holds, or if (4.2) (resp. (4.6)) holds. For this one has to compute the 
ranks of the principal submatriccs Ms{y) of Mt{y) for s < t. Checking the 
rank of a matrix consisting of numerical values is computationally sensitive. 
This is carried out using singular value decomposition which at the same 
time can be used to generate a basis of the column space; sec the next section 
for more details. The determination of the rank is done by detecting zero 
singular values or a decay of more than le-3 between two subsequent values, 
where singular values less than le-8 are declared to be zero. 

4.4.3 Computing a basis for the column space of Ms_i(y) and the 
quotient space K[a;]/J 

We indicate here how to compute a basis of the column space of the matrix 
Ms-i{y). Under some conditions (recall Corollary 3.8), such basis also yields 
a basis of the quotient space lK[x]/J (as before, J := (KerMs(y))), which 
is needed for the computation of the multiplication matrices. The choice of 
this basis will have an influence on the numerical stability of the extracted 
set W of solutions and on the properties of the basis for J as well. 

Using singular value decomposition. It is a well known fact from linear 
algebra that a numerically stable way of finding an orthonormal basis B for 
the column space of a matrix M is to use its singular value decomposition 
(SVD): M = UT,V*, where U, V arc unitary and S is diagonal with nonneg- 
ative entries. The diagonal entries of T, are the singular values of M (i.e., the 
square roots of the eigenvalues of MM*); the number of nonzero diagonal 
entries of S is thus equal to r := rank M. Then the set {Ui, . . . ,Ur} of 
columns of U corresponding to the nonzero diagonal entries of S forms an 
orthonormal basis of the column space of M. As we already did perform a 
SVD to determine the rank of the matrix M := Ms-i(y), this computation 
comes with no extra effort. For i = 1, ... ,r, lei bi := CT-i,x^i be the poly- 
nomial with vector of coefficients vec(6j) = Ui. The next lemma shows that, 
under some rank condition, {bi, . . . , 6,.} is a basis of K[x]/J. 

Lemma 4.16. Let {Ui, . . . , Ur} be a basis of the column space of Mg-iiy), 
let bi := Cj-i^x^i (i = 1, . . . ,r), and assume that the rank condition (4-1) or 
(4-5) holds. Then the set {bi, . . . , 6^} is a basis of the quotient space K.[x]/J. 

Proof. As the rank condition (4.1) or (4.5) holds, we know from Corol- 
lary 3.8 that dimIK[,T]/J = r := rankM^_i(y) and JnlKs_i[x] = KerM^_i(y). 
Hence it sufliccs to show that {bi, . . . , br} is linearly independent in ]K[x]/J. 
For this assume Ajbi G J, i.e., Cl-i,xi'^i=i ^i^i) ^ "^^^ vector 

p := Yll=i ^i^i li^s ™ the column space of Ms_i(y). On the other hand, 
p G KerMs_i(y) since the corresponding polynomial p lies in J fl 
Therefore, p = which implies that all Aj = 0. □ 

Using a 'greedy' algorithm. If we want to compute a border basis or a 
Grobner basis with our algorithm, we need a monomial basis B (i.e., B C T^) 
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for the quotient space ]K[a:;]/J. For this, it suffices to construct a set B 
indexing a maximum principal nonsingular submatrix of Ms_i(y) as B is 
then a basis of J under the conditions of Corollary 3.8. One can apply 
the following simple procedure (proposed in [29]) for constructing B: Scan 
monomials in T^^s-i by increasing degree, starting with to = Ij = a^i, ^2 = 
X2, ■ ■ ■■ Initialize B := {to}. Let B be the current set and tk be the current 
monomial to be scanned. If ;S U {tk} indexes a linearly independent set 
of columns of Ms_i(y), then reset B := B U {tk}, otherwise scan the next 
monomial tk+i ■ This procedure is 'greedy' in the sense that one keeps adding 
as many low degree monomials as possible to the basis. One can stop as 
soon as \B\ = r. Alternatively, one may construct a reduced row echelon 
form of Ms-i{y) using Gauss Jordan elimination with partial pivoting (pivot 
variables serve as basis B), see [21]. One can verify afterwards whether the 
constructed basis B is an order ideal; it turns out that this is the case in 
most tested instances. 

The greedy sieve algorithm described earlier in Section 2.4 produces 
directly an order ideal basis. Indeed, given any graded monomial ordering 
y, we can apply it to obtain the set B = By of standard monomials for this 
ordering, forming an order ideal basis of ]K[a;]/J (as we know from Lemma 
2.5 that By is contained in Tn,s-i under the conditions of Corollary 3.8.) See 
Section 4.4.5 for an extension to the case of an arbitrary monomial ordering. 

Note that, although desirable from an algebraic point of view, monomial 
bases for ]K[a;]/J sometimes lead to a less accurate set W of extracted so- 
lutions as compared to those extracted with a polynomial basis B based on 
SVD; see e.g. Examples 5.4, 5.5. 

4.4.4 Computing formal multiplication matrices 

Let B = {bi, . . . , bj.} be a basis of the column space of Ms_i(?/). By the as- 
sumptions of Proposition 4.1 or under the rank conditions (4.1) or (4.5), 
there exist scalars a^!"''' {k = l,...,r) for which Xibj — Ylk=i '^k''^ £ 
KerMs(y), for all i = l,...,n, j = l,...,r. Then the vector (a^*-'^)^^j^ 
is the jth column of the (formal) multiplication matrix Xxi- We indicate 
how to compute from Ms{y). 

Suppose first that ;S is a monomial basis, i.e., B C Tn,s-i- Let Mq denote 
the principal submatrix of Mg (y) indexed by B and let P^, . be the submatrix 
of Ms{y) whose rows are indexed by B and whose columns are indexed by 
the set XiB := {xibj \ j = 1, . . . ,r}. As observed in [29], we have 

X,^ = M^'P,^. (4.11) 
Indeed, for b G T^^s, let Cf, denote the column of Ms{y) indexed by b re- 
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stricted to the rows indexed by B. Then, 

r 

C.,bj=^a'i'^a,=MBa^''^, (4.12) 

k=l 

i.e., a(*^) = M-^C^^h,, which gives = M-^P^^. 

Suppose now that ;S is a polynomial basis obtained via SVD, as explained 
above. That is, bi = JJi where {C/i, . . . , Ur} is an orthonormal basis of 
the column space of Ms-\{y) and thus of Ms{y) under the rank condition 
(4.1) or (4.5). As in the monomial case, the formal multiplication matrices 
can be derived from Ms{y). Let Px^ denote the submatrix of Ms{y) with 
columns indexed by Xifn^s-i and with rows indexed by Tn,s-i- Let U denote 
the matrix with columns Ui, . . . ,Ur, and set Pj,. := U^Px^U and Mg := 

Ms-i{y)U . Then, Mjs is nonsingular. Moreover, 

MBXxi=Pxi, (4.13) 

which allows the computation of = where S is the 

diagonal matrix containing the positive singular values of Ms_i(y). We ver- 
ify that (4.13) holds. By construction, the polynomial Xibj — "^^.^i ci'^ ^'^k ~ 
Xi(^_i JJj - YJk^i a^f \'^_^^JJk lies in KerMs(y). This implies = Px^Uj - 
Ms-i{y)Ua^'^^ and thus U^Px^Uj = U'^ Ms-i{y)Ua^^^^ = MBa^'^\ which 
shows that the two matrices and M^^xi have identical jth columns. 

4.4.5 Constructing a basis for the ideal J := {Ker Ms{y)) 

A linear basis of KerM^ (y) . The simplest way of producing a basis for the 
ideal J = (KerMs(y)) is simply by considering a linear basis of KerMs(y). 
Such a basis can be found by using again a singular value decomposition 
for Ms{y). Indeed, if Ms{y) = U^V* is the SVD, then the columns Vi of 
V corresponding to the zero diagonal entries of S (the zero singular values 
of Ms{y)) form an orthonormal basis of KerMs(y). Then the polynomials 
CTx^i corresponding to the zero singular values of Ms{y) form a basis of J. 
A drawback of this basis however is that it is usually highly overdetermined 
and has a large cardinality, equal to IT^^sl — TajakMs{y). 

A border basis. As shown in [40, Sec. 8.2, Ch. 10], it is desirable to 
avoid overdetermined bases for J because it could lead to inconsistencies in 
the basis for numerical reasons. To avoid this drawback, border bases are 
proposed in [40] and their numerical properties are investigated. If during 
the construction of the formal multiplication matrices an order ideal basis 
B of IK [x]/ J was used, we deduce immediately a border basis consisting of 
the polynomials 

r 

Xibj - ^ al^hk for Xibj G dB. (4.14) 
k=l 
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A Grobner basis. If the monomial basis B of ]K[a;]/ J is the set of standard 
monomials By with respect to a monomial ordering y obtained, e.g., with 
the greedy sieve algorithm, then the border basis in (4.14) is actually a 
Grobner basis with respect to the monomial ordering When ;^ is a graded 
monomial ordering then, in view of Lemma 2.5, By. C T„ ^ and thus By can 
be found with Algorithm 1 applied to {I{Vk{I)),^, s), using the following 
independence oracle: a subset of T^^^ is independent in K[x]/I{Vk{I)) if and 
only if it indexes an independent set of columns of Ms{y). In general when 
>- is not a graded monomial ordering we are not assured to find By within 
Tn^s- However we can proceed as follows. As Ms{y) is a flat extension 
of Ms_i(y), by the results in Section 3.3, there exists an extension y G 
M"^" (for any t > s) such that Mt{y) is a fiat extension of Ms{y). As 
(KerMs(y)) = (KerMt(y)) = I{Ve.{I)), a subset of Tn,t is independent in 
K[x]//(ViK(-^)) if and only if it indexes an independent set of columns of 
Mt{y). Thus to find By, apply Algorithm 1 iteratively toi = s + l,s + 2,... 
until finding Bt = Bt+i, in which case we know from Lemma 2.6 (ii) that 
By = Bt- Remains only to address the question on how to find the fiat 
extension y. The existence proof in [12] is constructive and can roughly be 
sketched as follows (see also [29] for details). Say we want to construct a flat 
extension C := Ms+i{y) of B := Ms{y), under the assumption that B is a 
fiat extension of Ms-i{y). We indicate how to construct the column C(-,7) 
of C indexed by a monomial of degree s + 1. Write, say, x'^ = xix^. By 
assumption, the column B{-,(5) of B indexed by x^^ can be expressed as a 
linear combination X]jQ,|<3_i Aa5(-, a) of columns indexed by T„,s-i; then 
define C(-,7) as X^|q,|<s_i \aC{-,a + e^). Note that this construction relies 
on the fact that the kernel of moment matrices enjoys ideal-like properties. 

4.4.6 Summary of the algorithm 

Algorithm 2 below summarizes our algorithm. This algorithm has been im- 
plemented in Matlab using the Yalmip toolbox [30]. For solving the semidef- 
inite program (4.9) the semidefinite solver SeDuMi-1.05 [41, 42] is used. As 
described above and can be seen in the examples in the next section, the 
rank detection is the most critical task. This was the main motivation for 
the weaker conditions from Section 4.1, which extend the possibility of ex- 
tracting solutions. In the examples below this deficiency is clearly indicated 
in some rank sequences not exactly matching the theory. 

5 Numerical Examples 

We present here the results of our algorithm applied to some examples, 
mainly taken from the literature. In each example, we specify the ideal / 
by its generators hi,. . . ,hm- Let us explain Tables 1-5 below. At a given 
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Algorithm 2 Numerical border basis computation: 

Input: Polynomial generators hi for I := {hi, . . . hm) ^ K[x] and relaxation 

order t G N 

Output: A basis for an ideal J C I(yK{I)), the set Vc{J), and a basis B 
for the quotient ring K[x]/J 

1: Solve the SDP (4.9). If the SDP is infeasible, return Vk{I) = 0. Other- 
wise, return a feasible solution y for which Mt{y) has maximum rank 

2: Compute SVD for all principal submatriccs Ms{y) (s = 1, . . . , t) 

3: Determine rankMs(y) (s = 1, . . . ,i) and check whether the conditions 
of Prop. 4.1-4.11 hold 

4: Fix s (for which one of Prop. 4.1 - 4.11 applies) 

5: Compute a basis B of the column space of Ms-i{y): 

a) using the SVD decomposition {B is a polynomial basis) 

b) using a greedy algorithm {B is a monomial basis) 

c) using a greedy sieve algorithm {B is the set of standard monomials 
for a monomial ordering >-) 

6: Compute the multiplication matrices = ^s^Pxi 
7: Compute a basis for the ideal J 

a) a SVD basis of Kci Mgiy) (requires rankM5(?/) = rankM5_i(y)) 

b) a border basis of ./ (requires that S is a monomial basis) 

c) a Grobner basis (requires that B is the corresponding set of standard 
monomials) 

8: if the conditions of Prop. 4.1 are met then 

9: return a border /Grobner basis of J C I{y^[I)), a basis B of ]K[a;]/J, 

and the set Vc{J) 
10: else 

11: return ERROR: "No extraction possible. Increase relaxation order 
t." 

12: end if 

Remark 4.17. In view of Propositions 4.6, 4.12, the algorithm terminates 
for t large enough and finds J = /(Vr(/)) or I{yc{I)). The algorithm can 
also be used for testing existence of solutions. Let us give some details e.g. 
in the real case. If at step 1 one detects infeasibility of the SDP then one can 
already conclude Vr{I) = 0. Suppose now the SDP is feasible. At step 4, as 
observed in Remark 4.8, the conditions of Proposition 4.4 cannot be met if 
Vw{I) = 0, but it could be that the conditions of Proposition 4.1 or 4.2 are 
met. In that case one can extract a set W D Vm.{I). Then one can simply 
test whether the points of W satisfy the given equations hi = . . . = hm = 
to detect whether Mr(/) is empty or not. 
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order t, let y be the optimal solution to (4.9) returned by the SDP solver. 
The abbreviations 'MON' and 'SVD' refer to using a monomial base of the 
quotient space, or a base found via the SVD method. 

• The column 'rank sequence' shows (rankMo(y), . . . , rankMj(y)). 

• The column 'extract, order' shows some numbers 'Smon('"mon)/ssvd(''svd)- 
When using a monomial base, Tmon is the smallest order at which the ex- 
traction procedure could be carried out and Smon is the order at which it 
was effectively carried out and gave the results reported here; analogously 
with the SVD method. 

• The column 'accuracy' shows the accuracy of the returned solutions, i.e., 
maxj^x \hj{x)\, where hj runs over the generators of I and x over the ex- 
tracted solutions. 

• The column 'comm. error' shows the commutativity error for the mul- 
tiplication matrices, i.e., max"^-^]^ abs(Xc. X^^. — Xx^X^^) (where abs(M) is 
the maximum absolute value of the entries of a matrix M). If the param- 
eter 'comm. error' is more than le-2, the multiplication matrices do not 
commute sufficiently and we then do not extract solutions. 

Example 5.1 This simple example from [11, p. 40] has two roots, both 
real. 

h2 = xlx2 - 2x\ 

/i3 = 2x|xi — x\ — 2X2 + ^1 



order 
t 


rank sequence 


extract, order 
MON/SVD 


accuracy 
MON/SVD 


comm. error 
MON/SVD 


3 


13 5 9 








4 


1 2 2 2 7 


4(2)/3(2) 


1.9717e-9/0.00013144 


9.676e-10/3.3908e-6 


5 


12 2 2 2 8 


4(2)/4(2) 


2.9557e-8/3.5325e-5 


1.8781e-ll/1.2291e-6 



Table 1: Results for Example 5.1 



Monomial basis of R[x]/I{Vr{I)): 

B = {l,xi}. 

Border basis for I(Vk(/)) (showing in bold the monomials in dB): 

91 = -X1+-X.I, 
g2 = -2x1 +X2, 
93 = -2x1 +X1X2. 
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Extracted real solutions Vm(/): 

xi = (2.12c-8,1.91e-6), 
X2 = (1,2). 

The first two polynomials gi , g2 of the extracted border basis form a reduced 
Grobner basis with respect to the graded reversed term order with xi -< X2- 
The basis of ^/l (= I{Vr{I)) as all roots are real) given in [11] has the form: 

2x2X1 - xl- 2x2 + ^1' xii^i - 1), a^2(-2 + X2)} 

and is obtained via Seidenberg's method described in the paragraph 'Related 
literature' in the Introduction. Computing a graded reversed Grobner Basis 
of vT (using tdeg in Maple) leads again to the set {51,(72} found by our 
method. Thus our method finds here a simpler set of generators for s/l than 
the classical method of Seidenberg. 

Example 5.2 This example is taken from the polynomial testsuite [6] (see 
http : //www-sop . inria . f r/saga/POL/BASE/2 . multipol/bif urc . html). It 

has 20 complex solutions among which 8 are real. This example illustrates 
the possibility of extracting solutions based on Proposition 4.1 in case none 
of the rank conditions are satisfied. 

hi = 5xi — 6x^X2 + X1X2 + 2x1X3 
h2 = -2x\x2 + 2x1x1 + 2x2x3 
h3 = x\ + xl- 0.265625 



order 
t 


rank sequence 


extract, order 
MON/SVD 


accuracy 
MON/SVD 


comm. error 
MON/SVD 


5 


1 4 8 16 25 34 








6 


1 3 9 15 22 26 32 








7 


1 3 8 10 12 16 20 24 


3(3)/-(-) 


0.12786/— 


0.00019754/— 


8 


1 4 8 8 8 12 16 20 24 


4(3)/3(3) 


4.6789e-5/0.00013406 


4.7073e-5/0.00075005 



Table 2: Results for Example 5.2 
Monomial basis of W^x]/ I{V^{I)): 



B = {1, Xi,X2, X3,xl, X1X2, X1X3, X2X3} . 
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Border basis of 



91 = 


— U.zo47yxi + U.44iz4a;ia;2 — i.oAUdxiXs + , 


92 = 


-1. 7276x3 - 0.080949x1 + 8.1433x2X3 + XiX2 , 


93 = 


-0.28763x3 - 0.0010314X? + 0.48126x2X3 + xg , 


94. = 


-0.0015073x1 + 0.01299x1X2 - 0.12111xiX3 + xixaxg 


95 = 


-0.26563 + X? + x| , 


96 = 


0.019164x1 - 0.44124x1X2 + 1.5403xiX3 + xix| , 


97 = 


0.022008x3 + 0.0010314X? - 0.48126x2X3 + x^xg , 


98 = 


-0.0018637x3 - 0.00067043X? + 0.026066x2X3 + x§ , 


99 = 


-0.00015166x1 + 0.00025958x1X3 + xix§ , 


910 = 


-0.0017335x3 + 0.01615x2X3 + X2X§ . 



Extracted real solutions: 

xi = (-0.515,-0.000153,-0.0124) , 

X2 = (-0.502,0.119,0.0124), 

X3 = (0.502,0.119,0.0124) , 

X4 = (0.515, -0.000185, -0.0125) , 

X5 = (0.262,0.444,-0.0132) , 

X6 = (-2.07e-5, 0.515, -1.27e-6), 

X7 = (-0.262, 0.444, -0.0132) , 

X8 = (-1.05e-5, -0.515, -7.56e-7) . 

Example 5.3 We now give an example for finding I{Vm.{I) fl S), where 
/ = {hi, . . . , h^) with 

/ll = Xl + X2 — 2 , 
h2 = X1X3 + X2X4 , 

/13 = Xix| + X2X4 - ^ , 
h4 = Xix| + X2X4 

and S is defined by the polynomial inequalities —1 < xi, X2, X3, X4 < 1. This 
example, taken from [45] , represents a Gaussian quadrature formula with 2 
weights and 2 knots, where one is interested only in the roots lying in the 

box [-1,+1]. 

Monomial basis of M[x]//(Vir(/) n S): 

B = {l,xs}. 
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order t 


rank sequence 


extract, order 
MON/SVD 


accuracy 
MON/SVD 


comm. error 
MON/SVD 


2 


1 2 11 


2(2)/-(-) 


0.00010224/— 


1.1124C-9/— 


3 


1 2 2 18 


2(2)/2(2) 


1.8985e-14/5.1015e-14 


1.2212e-15/1.4155e-15 


4 


1 2 2 2 24 


2(2)/2(2) 


3.5527e-15/8.5487e-15 


2.2204e-16/1.1102e-16 



Table 3: Results for Example 5.3 



Border basis of I{Vm{I) n S): 



91 


= -1 + xi 


92 


= -X3 + X1X3 


93 


= -1 + X2 


94 


= -X3 + X2X3 


95 


= -0.33333 + x| 


96 


= ^3 + X4 


97 


= 0.33333 + X3X4 



For this example the border basis is in fact a Grobner basis, e.g. with respect 
to a graded lexicographic order with xi )^ X2 )^ X4 )^ X3. Note however that 
it is not a reduced Grobner basis since 52 , 54 and 57 are redundant for this 
ordering. 

Extracted real solutions Vr{I) Ci S: 

xi = (1,1,-0.577,0.577), 
X2 = (1,1,0.577,-0.577). 

Example 5.4 This example: Katsura 5, is an example in with 32 com- 
plex roots, including 12 real roots. It is taken from 

http : / /www . mat . univie . ac . at/~neum/ glopt/ coconut /Benchmark/LibraryS/katsuraS . mod. 

^1 = 2x6 + 2a:i + 2x1 + 2x1 + 2x1 + xj - xi , 

/l2 = XgXs + X5X4 + 2X4X3 + 2X3X2 + 2X2X1 — X2 , 
/l3 = 2X6X4 + 2X5X3 + 2X4X2 + X2 + 2X3X1 — X3 , 
/14 = 2X6X3 + 2X5X2 + 2X3X2 + 2X4X1 — X4 , 
h5 = x'^ + 2xqXi + 2x5X1 + 2x4X1 — X5 , 
he = 2x6 + 2x5 + 2x4 + 2x3 + 2x2 + xi - 1. 

We cannot extract solutions with a monomial base since the multiplication 
matrices do not commute, but we can extract the following real solutions 
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order 
t 


rank sequence 


extract, order 
MON/SVD 


accuracy 
MON/SVD 


comm. error 
MON/SVD 


1 


1 7 








2 


1 6 16 








3 


1 6 12 12 


-/3(3) 


— /1.1928e-005 


— /2.3073e-007 



Table 4: Results for Example 5.4 



using a SVD basis: 



Xl 


= (0.277, 0.226, 0.162, 0.0858, 0.0115, -0.124) , 


X2 


= (0.59, 0.0422, 0.327, -0.0642, -0.0874, -0.0132) , 


X3 


= (1, -2.8C-7, 4.7e-7, 8.81e-7, -2.79c-6, -3.69e-6) , 




= (0.239, 0.0608, -0.0622, -0.0233, 0.186, 0.219) , 


X5 


= (0.441, 0.151, 0.0225, 0.219, 0.0935, -0.207) , 


Xq 


= (0.726, -0.0503, 0.122, 0.164, 0.11, -0.208) , 


XT 


= (0.462, 0.309, 0.0553, -0.102, -0.0844, 0.0917) , 


X8 


= (0.292, -0.101, 0.181, -0.0591, 0.193, 0.141) , 


Xg 


= (0.753, 0.0532, 0.191, -0.114, -0.146, 0.139) , 


xw 


= (0.409, -0.0732, 0.0657, -0.127, 0.252, 0.178) , 


X\l 


= (0.68, 0.266, -0.154, 0.0323, 0.0897, -0.0735) , 


Xl2 


= (0.136, 0.0428, 0.0417, 0.0404, 0.0964, 0.211) , 


Example 5.5 


The following example shows the computation of the radical 



ideal using complex moment matrices. 



hi = xl + X2 + X3 + l, 
h2 = Xl + xl + Xs + 1 , 
hs = Xl + X2 + xl + 1 . 

This ideal is not radical and admits 7 solutions, among which the solution 
(-1,-1,-1) has multiplicity two. We solve the SDP program based on the 
set Kf'^ (thus using full complex moment matrices) . The rank sequence for 
full and pruned matrices Mf{y) and Mf~'{y) are shown in Table 5. 
Monomial basis of C[x]//(Vc(/)): 

B = {1, Xi,X2, X3, XiX2, X1X3, X2X3}. 
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order 


rank sequence 


extract, order 


accuracy 


comm. error 


t 


Mf (y), (M2€(y)) 


MON/SVD 


MON/SVD 


MON/SVD 


1 


(1 4), 










(1 7) 








2 


(1 4 7), 










(17 7) 








3 


(1 4 7 7), 


3(3)/3(3) 


0.0005719/0.00022538 


0.00041241/0.00043871 




(1 7 7 7) 









Table 5: Results for Example 5.5 



Border basis of /(Vc(-^)): 
51 = 1 + a;2 + X3 + Xi 

g2 = -1- Xi+ X2- X3 + X2Xs + X^X2 , 

53 = -1 - Xi - X2 + X3 + X2X3 + X1X3 , 

54 = 3.9993 - 0.99984x1X2 - 0.99984xiX3 - 0.99984x2X3 + X1X2X3 , 

55 = 1 + Xl + X3 + x| , 

g6 = -I + X1-X2- xs + X1X3 + xix| , 

57 = -1 - Xi - X2 + X3 + X1X3 + x|x3 , 

58 = 1 + a;i + X2 + x§ , 

59 = -1 + Xi - X2 - 0^3 + X1X2 + Xix| , 
510 = -1 - Xi + X2 - X3 + X1X2 + X2X§ . 

Extracted solutions (via the monomial basis): 

xi = (-1 - 8.15e-lli, -1 + 4.37e-lli, -1 - 4.24e-12i) « (-1, -1, -1) , 
X2 = (-1.16e-5 + 1.4H, 0.999 - 1.4H, 0.000654 + 1.4H) « {V2i, 1 - V2i, V2i) , 
X3 = (-4.8e-5 - 1.4H, 1 + 1.4H, 0.000147 - 1.4H) « (-\/2i, 1 + V2i, -V2i) , 
X4 = (-9.92e-7 + 1.4H, 0.000713 + 1.41i, 0.999 - 1.41i) « {V2i, V2i, 1 - V2i) , 
X5 = (-3.98e-5 - lAli, 0.000146 - 1.41i, 1 + 1.41i) « {-V2i, -V2i, 1 + V2i) , 
X6 = (1 + 1.4H, -0.000149 - 1.4H, -0.000145 - 1.4H) ^ (1 + V2i, -V2i, -V2i) , 
X7 = (1 - 1.41i, 0.000103 + 1.41i, 0.000104 + 1.41i) « (1 - V2i, V2i, V2i) . 



40 



For this example more accurate solutions: 

xi = (-1 - 2.51e-llz, -1 - 9.85e-llz, -1 - 5.99e-llz) , 
X2 = (-1.73e-5 + l.m, 1 - 1.4H, -1.29e-5 + lAli) , 
= (-3.78e-5 - lAli, 1 + 1.41i, -8.52e-5 - 1.41z) . 

X4 = (-2.83e-5 + 1.41f, 4.23e-5 + 1.41i, 1 - 1.41i) , 
X5 = (-4e-5 - 1.41i, 3.88e-5 - 1.41i, 1 + 1.4H) , 
X6 = (1 + 1.41i, -0.000191 - 1.4H, 6.2e-5 - 1.4H) , 
X7 = {1- 1.41i, -4.61C-5 + 1.4H, 0.000104 + 1.41i) , 

could be obtain by means of the SVD-method. 

6 Concluding Remarks 

In this paper wc have provided a new semidefinite characterization of the 
real radical ideal of an ideal / C M.[x] as well as an algorithm to compute all 
(finitely many) points of and a set of generators (or a Grobner base) 

of /(T^(/)). The main feature of our approach is its real algebraic nature 
as it avoids considering complex zeros, and does not need to compute a 
Grobner base of /. An essential step in our algorithm consists in solving the 
semidefinite program (4.9). Thus our algorithm is numerical. 

Let us briefly mention the numerical versus numeric-symbolic (or arbi- 
trary precision) issue. Some advocate that only computation with arbitrary 
or guaranteed precision should be permitted while others admit some nu- 
merical imprecision; see e.g. Revol and Rouillier [35], Stetter [40]. Clearly, 
being numerical in nature, the algorithm of the present paper admits some 
intrinsic numerical imprecision, no matter how good are (or will be) the 
SDP software packages. At this stage, the only answer we propose in this 
'approximate vs exact' debate is to validate or invalidate the method by ex- 
periments. For instance, on a significant sample, compute J 7(1^(7)) with 
our method and check afterwards by symbolic methods if Vc{J) = V^{I). 
On the other hand, the present algorithm is rather intended to illustrate 
that the new semidefinite characterizations of Vk(/) and I(V%(I)) are di- 
rectly implementable in a relatively simple manner; clearly, its numerical 
features (like precision and stability) need further investigations beyond the 
scope of the present paper. 
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